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Abstract 

We analyze the spectrum of excitations around the Gubser-Klebanov-Polyakov (GKP) rotating 
string in the long string limit and construct a parametric representation for their dispersion relations 
at any value of the string tension. On the gauge theory side of the AdS/CFT correspondence, i.e., 
in the planar = 4 Super- Yang-Mills theory, the problem is equivalent to finding the spectrum of 
scaling dimensions of large spin, single-trace operators. Their scaling dimensions are obtained from 
the analysis of the Beisert-Staudacher asymptotic Bethe ansatz equations, which are believed to solve 
the spectral problem of the planar gauge theory. We examine the resulting dispersion relations in 
various kinematical regimes, both at weak and strong coupling, and detail the matching with the 
Frolov-Tseytlin spectrum of transverse fluctuations of the long GKP string. At a more dynamical 
level, we identify the mechanism for the restoration of the SO (6) symmetry, initially broken by 
the choice of the Berenstein-Maldacena-Nastase vacuum in the Bethe ansatz solution to the mixing 
problem. 
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1 Introduction 



The AdS/CFT correspondence [Ij has considerably deepened our understanding of the dynamics of 
strongly coupled gauge theories and shed light on the longstanding problem of constucting their string 
theory duals. In the most well studied example of the maximally supersymmetric extension of Yang- 
Mills theory (aka N = A SYM theory), that we are going to consider here, it offers a geometrical 
interpretation of the gauge dynamics in terms of strings propagating on AdSs x background. One 
of the illustrations of how this dictionary works is well established for the spectrum of scaling dimen- 
sions of composite operators in the conformal SYM theory. In the planar limit, the scaling dimensions 
are complicated functions of the 't Hooft coupling constant = A/lGvr^, with A = g'^^^Nc, and can be 
computed at weak coupling order by order in the perturbative expansion in g^. The prediction from 
the AdS/CFT correspondence is that this spectrum of scaling dimensions simultaneously describes 
the spectrum of energies of a single string propagating freely on the background geometry. This in- 
terpretation becomes even more manifest at strong coupling, i.e., at large string tension g ^ 1, where 
the quantum fluctuations of the string are suppressed. This conjecture has triggered an enormous 
flood of theoretical activity over the past decade and has led to dramatic improvement of our under- 
standing of theories on both sides of the duality. For the spectrum of scaling dimensions/energies, it 
has culminated with a proposal of a set of equations [21 [3l lU [5] which implements hidden integrable 
structures of the gauge/string system [SI El El El UHl E] and which is believed to solve exactly the 
spectral problem for any value of the coupling constant g. In this paper, we will apply these equations 
to unravel the spectrum of excitations around the Gubser-Klebanov-Polyakov (GKP) long rotating 
string p!2] and construct their dispersion relations at any value of the coupling constant. 

There is already a wealth of insights both at weak and strong coupling for the problem that will be 
the focus of this work. On the string theory side, the perturbative spectrum of fluctuations around 
the long GKP string, as well as their leading-order dispersion relations, have been known for a long 
time [131 E] • It is also known that this spectrum gets corrected at the non-perturbative level with 
the massless modes of the string acquiring a mass by dimensional transmutation ^14j . At the classical 
level, but now in the regime of giant excitations, the dispersion relation was recently constructed in |15j 
using the finite gap machinery and an interesting aspect regarding its large-momentum asymptotics 
was notably uncovered. On the gauge theory side of the duality, the problem maps into the analysis 
of the spectrum of scaling dimensions of large spin operators, for which many features have been 
already elucidated [TUl |T71 [TH [T21 120] and for which many results were already extended to all 
loops. In this paper, we would like to complete the picture and present a comprehensive account of 
the spectrum of excitations and of their dispersion relations. 

This investigation is not purely academic and was greatly motivated by the recent analysis of [21j 
where it is shown to be directly relevant to the computation of on-shell scattering amplitudes. The 
bridge between the two comes from their respective relationship with the cusped light-like Wilson 
loops, which has been known for some time for the large spin spectral problem [221 flSl [23l and 
uncovered more recently for scattering amplitudes [25]. This gives hope that the considerations 
and results deduced in this paper could be helpful to exploring the relation between the integrable 
structures observed in the study of scattering amplitudes at strong coupling |26j and the ones that 
govern the spectrum of scaling dimensions. 

Getting back to the problem at hand, it is worth recalling some of the remarkable features of 
the GKP string. By construction it describes the leading Regge trajectory of string eigenstates 
carrying single 50(2,4) spin S and minimal global-time energy A, and it can be thought as the 
AdSs analogue of the Berenstein-Maldacena-Nastase (BMN) single SO (6) spin J solution [271 112j . 
which is the standard reference state for building up the spectrum of scaling dimensions / energies using 
integrability [8l|9l[l0]. Semiclassically {g ^> 1) the GKP string is described by a folded string rotating 
around its center of mass in AdSs C AdSs x S^. At generic values of the spin, it defines a complicated 
solution to the classical string equations and thus represents itself an intricate background for the 
semiclassical expansion of the string cj-model [13J. In the large spin limit however, the string gets 
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long, with proper length ~ 2 log 5111 ultimately stretching all the way to the boundary of space at 
S = CO, and becomes homogeneous. Its proper energjH ^vacuum is then uniformely distributed and 
scales with the length leading to the logarithmic scaling [12] 

^vacuum = A - 5 = ^2(5) log S + 0{log^ S) , (1.1) 

where A2{g)/2 stands for the energy per unit of length given classically by the string tension A2{g)/2 = 
2g + . . . = \/A/27r + . . . . The scaling p.l|) persists at higher loops and its string 1/g expansion 
is known explicitely up to two loops \13\ I28j . A natural question arises what are the excitations 
propagating in this background and what are their dispersion relations. As we already alluded above, 
the answer is known to leading order at strong coupling and can be obtained by considering the 
spectrum of transverse fluctuations of the string [13i, il4j . It was found that the dispersion relations 
are all relativities and that the elementary world-sheet excitations are 5 massless bosons for rotations 
in S^, 2 mass-\/2 bosons for rotations of AdSs in AdSs, 1 mass-2 boson for the transverse fluctuation 
in AdSs, and finally 8 mass-1 fermions. The energy of the state made out of M excitations is then 
given bjo 

M 

E = A-S = ^vacuum + Ej{pj) + . . . , (1.2) 

i=i 

with pj = TTTLj/ log S for some mode numbers rij, and where the ellipsis stands for l/S'-suppressed 
corrections. It remains a challenge to compute corrections to this spectrum by conventional c-model 
methods. The AdS/CFT correspondence hints that we can solve this problem by considering the 
spectrum of scaling dimensions of large spin operators. Integrability suggests that the decomposi- 
tion ()1.2p should persist, in the interacting theory, giving the means to extract individual energies of 
the excitations. 

The gauge theory dual of the GKP string is a single-trace operator belonging to the (universal) 
twist-two multiplet [12j , though it is not entirely clear (to the author at least) what precise gauge- 
theory operator it corresponds to|l| For our purposes, it is not relevant which twist-two operator maps 
exactly to GKP string, and in practice it is more convenient to work with the following superconformal 
descendant!! 

vacuum = GKP string = tr ZDIz + ... , (1.3) 

carrying the Lorentz spin S, and where Z is one of the three complex scalars of the gauge theory, 
D+ = n'^Df, = Do + Ds is the light -cone covariant derivative, and dots indicate that the form of 
the operator is renormalized. The operator representation (II. Sp obscures however the fact that the 
GKP vacuum is a SO{6) singlet. An immediate consequence of the mapping p.3p is the recovering 
of the logarithmic scaling (jl.ip to all loops [22], with A in (jl.ip the scaling dimension and with 
^2(5) = 2rcusp(5') = 4(7^ -|- 0{g^) being twice the cusp anomalous dimension [SSUSO]. Since 2 log 5 
stands for the physical length in the game|£| the identification (jl.ip implies that the cusp anomalous 
dimension has the meaning of the vacuum energy density |14j . 

The next step is to add excitations on top of the ground-state operator (jl.Sp and extract their 
energies using the decomposition (jl.2p . A natural way of doing this is to allow for operator insertions 

^With the radius of curvature of the warped geometry normalized to one. 

^That is the energy measured in the rotating frame of the long GKP string. 

^Perturbative excitations are free to leading order at strong coupling. 

possible candidate is the superconformal primary which is a 50(6) singlet. 

^The spin of the operator may differ from the spin of the GKP string since we are considering a descendant. This 
leads to a finite shift that is irrelevant at large spin within the accuracy of our approximation. 

®The interpretation of the quantity 2 logs' as defining the physical length in the problem we are considering is 
not straightforward at weak coupling. It is however possible to construct a geometrical picture for the large spin 
operator (|1.3p in which this identification becomes transparent [14|. This picture is obtained by conformal mapping 
of the gauge theory from R'^''^ to AdSs x S^. The twist-two operator then translates into a flux conflguration which is 
stretched along a particular direction with physical size 2 log S. Exciting this conflguration is equivalent to exciting the 
operator (|1.3p and dual to exciting the GKP string. 
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m (mill For a single excitation, for instance, one expects the following identity to take place, 

one-particle state = tr ZDI^'^dI^'Z + ... , (1.4) 

where <I> stands for a generic local operator of the gauge theory and with S ^ Si + S2 ^ 1. Under 
renormalization, the operator p.4|) can mix with similar operators that differ in their relative numbers 
iSi — 52 of covariant derivatives. The interpretation is that the operator $ propagates through the 
background of covariant derivatives. It carries a momentum p which, roughly, accounts for the 
inequivalent ways of distributing the covariant derivatives around it. Imposing the condition that 
the operator (|1.4p carries a definite scaling dimension A should fix its form and translates into a 
quantization condition for the momentum p of <I>. 

Obviously, this picture is directly inspired from the treatment of excitations on top of the BMN 
vacuum tr Z"^ [27j . The difference here is that since the vacuum operator ()1.3p has a complicated 
structure, with an intricate mixing, it would be difficult to make use of it outside of the realm of 
integrability. Keeping this in mind, we can push the analogy with the BMN analysis even further 
and ask what are the (lightest) elementary excitations. For the BMN case, they are those that make 
the minimal contribution to the BMN energy A — J. For the GKP string, minimizing the energy 
E = A — S = twist + 0{g^) suggests that the elementary excitations are twist-one partons of the 
gauge theory. They are also known as light-cone operators [31] and are the building blocks of the 
quasi-partonic operators [32j . The pattern of elementary excitations is then composed of eight bosonic 
and eight fermionic operators that form multiplets under the residual 50(6) symmetry of the GKP 
vacuum: we get 6 (real) scalar fields in the 6 of su(4) = so(6), 4/4 twist-one components of left /right 
Weyl spinors in the 4/4, and 2 twist-one components of gluon field strength tensor in the Keeping 
only the su(4) highest-weight fields, the twist-one partons are 

F+x = «F^,, ^'+, Z, §+, = n'^nlF^, , (1.5) 

with n_L = (0, 1, i, 0), n_L = (0, 1,— i,0), the vectors transverse to the lightcone direction, n+ • n_L = 
n+ • n_\_ = 0, and with the -|- subscripts on the spinor fields indicating their twist-one component. 

At the tree level, the twist-one excitations form a supersymmetric pattern and their energy is 
given by their mass which is equal to their twist. Going beyond this order entails solving the mixing 
problem for the light-cone operators p.4|) . At one loop, it is equivalent to the diagonalization of the 
Hamiltonian for the s[(2) XXX Heisenberg spin chain, whose spin module depends on the flavor of 
the parton (see [33.j for the case at hand). It is analyzed systematically by means of a set of Bethe 
ansatz equations whose solution at large spin 5 has been given in [T7]. It is immediately found that 
the energy of a twist-one excitation is [17j □ 

E{u) = l + 2g'^ {tp {s + iu)+il){s- iu) - 2V'(1)) + 0{g^) , (1.6) 

where it is a rapidity that parameterizes the motion of the excitation through the "bath" of covariant 
derivatives, and where ip{z) = V'{z)/T{z) is the logarithmic derivative of the Euler Gamma function 
T{z). Here s stands for the (light-cone) conformal spin of the excitation and accounts for different 
flavors, with 2s = 1, 2, 3, for scalars, fermions and gauge fields, respectively. The analysis of [T^ also 
contains the effective Bethe ansatz equation that fix the rapidity of the excitation (once confined on 
the GKP string). From their expression, we read that the momentum p, conjugate to the GKP length 
2 log 5, is 

p{u) = 2u + 0{g^), (1.7) 

''it is necessary to increase the length of the operator if we want to put excitation on top of it. This is 

because all length-two operators fall into the same twist-two multiplet [341 [TO] and since taking descendant is a global 
transformation which 'does not propagate'. 

*The twist-one component is the one that carries maximal Lorentz spin S in a given S'0(3, 1) multiplet. Other 
components possess higher twist. 

^The result of [T7] differs actually from (|1.6p by a vacuum energy shift (namely the energy of the spin chain when no 
covariant derivatives act on it). This constant depends on the theory in which the conformal spin chain is embedded. 
For A/" = 4 SYM theory, see [lb] for instance, we need to add 4g^(i/)(2s) - ■0(1)) to the result of [TTj to obtain (fTB)) . 
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for any excitation!^ Combining (II. 6p and ()1.7p gives the sought one-loop dispersion relations. 

In this paper, we will compute subleading corrections to (|1.6U1.7p . and derive the all-loop dispersion 
relations from the Beisert-Staudacher asymptotic Bethe ansatz equations [2] . The representation for 
the dispersion relations is parametric and is expressed in terms of particular functions that solve the 
Beisert-Eden-Staudacher (BES) equation ^ for the ground-state distribution of roots. Our analysis 
follows the seminal papers [351 [T7\ \TE[ [TO] for dealing with the large spin limit of the Bethe ansatz 
equations, and, as already mentioned above, it overlaps with previous analyses. In the particular case 
of a scalar, for instance, subleading corrections to the energy (jl.6p were already computed in [36j 
up to three loops, by means of the asymptotic Baxter equation for the s[(2) sector [37j. It is also 
possible to derive the all-loop expression for the energy, given in this paper, directly from the analysis 
performed in [3H] of the Freyhult-Rej-Staudacher (FRS) equation [TU] for the so-called generalized 
scaling function [171 119j . Indeed, this function computes the ground-state energy density of a gas 
of scalars (confined on the GKP string) at finite density [13] and thus encodes the all-loop energy 
of a scalar. It has to be complemented by a similar expression for the momentum to get the full 
parameterization of the dispersion relation. 

As a particular application of our formulae, we can observe the merging of the spectrum of ele- 
mentary excitations of the gauge and string theories, see Fig. [TJ The fate of the scalars has been 
already elucidated in [13] • Their mass decreases with the coupling and becomes exponentially small 
at strong coupling. The perturbative string analysis detects only 5 massless excitations, which are 
Goldstone bosons for rotations in S^. But these are not the genuine asymptotic excitations and the 
actual spectrum contains 6 massive scalars with mass m ~ e~'^^, in agreement with the gauge theory 
expectations [13] • These features were already tested in the framework of integrability [Ml ESI SO] in 
the context of the FRS equation, and, not surprisingly, are also encoded in our relations. The case 
of fermions was also covered by the analysis of [13] in which it is explained that they ought to be 
exactly mass 1 excitations, as a consequence of the remnant supersymmetries spontaneously broken 
by the GKP vacuum. Here this property is also a consequence of supersymmetry and is tied to the 
multiplet joining/splitting [Si [2]. For the gauge fields, the weak coupling expansion indicates that 
their mass increases with the coupling and analysis at strong coupling shows that it settles down 
at \/2!l3 They find their place in the spectum of perturbative string fluctuations where they are 
identified with the two transverse mass \/2 bosons. For illustration of these considerations, we report 
the one-loop corrections to the masses of twist-one excitations, that immediately follow from (II. 6[ 
[LTP evaluated at u = o!ll 

^scalar = 1 " 8c/^ log 2 + 0(5^) < 1 , 

'^fermion — 1 ; (-^'8) 
mgauge field = 1 + 85^(1 - log 2) + 0(/) > 1 . 

It is not much difficult to compute corrections to the relativistic laws for the perturbative excitations 
at strong coupling, which read up to 0{l/g), 



with the masses 



(1.9) 



m-scalar = 0(1/5°°) , mfermion = 1 , "Igaugc field = ^2 - ^J^^ + 0(l/g^) . (1.10) 



^°Note that, dispite its appearance, the expression p.7p for the momentum is a genuine one-loop computation. 

^^Our analysis for the mass of gauge field overlaps with the one of |20) . in which the embedding of the gauge-field 
excitation in the aU-loop Bethe ansatz was uncovered (see also [41] for previous discussion of the procedure used in |20)1. 
Though the aU-loop integral representation is the same here and in [5^, the strong couphng expansions disagree. 

^^Note that all 8 fermions have the same mass, and idem for the 2 gauge fields, due to the 0(2) symmetry that rotates 
the _L directions and also exchanges left and right. 
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Figure 1: Spectrum of masses for gauge and string elementary excitations. The question mark is for 
the missing mass 2 boson predicted by the perturbative analysis in the string theory [T31 Its 
mass should not receive corrections |14J . There are candidates for it but their stability at non-zero 
momentum is not guaranteed, since decay channels into two fermions are allowed. 



The numerical coefficients entering (|1.9p at 0{l/g) are found to be 



Cscalar — 8(12)'^/^7r^ ' '^fermion — ^ i Cgauge-field — g • 

Looking further at the table in FigH] one observes that the twist-one excitations do not exhaust 
the perturbative spectrum of string fluctuations. The (heavy) mass-2 boson of the string theory is 
apparently missing. According to [H] its mass is protected and thus should be exactly 2. Twist- two 

candidates with correct quantum numbers exist for it, like = n'^n'^F^^, with n_ = (1,0,0, — 1), 

or possibly -D_ = ri- ■ D. These fields do not find room in our analysis as stable (all-loop) asymptotic 
excitations. The argument of [lA^ do not seem actually to guarantee that the boson is stable at non- 
zero momentum (except at strong coupling due to the boost invariance) and a possible decay into two 
fermions appears quite realistic, as already advocated in [21]. A similar phenomenon was actually 
observed in (42j in a relative set upEl Whatever the fate of the mass-2 boson can be, its presence 
raises the question of possible higher-twist excitations. There exist already evidences for them. They 
come from the spectroscopy of length-three operators [l3] that shows that higher-conformal-spin s[(2) 
spin chains are embedded into the super spin chain [S] of the planar gauge theory. The conformal spin 
s is then allowed to take arbitrarily large values 2s = 1, 2, 3, ... , the first three ones being realized by 
the twist-one fields (|1.5|) and the others by the (length-one) partons 

D'C'F+^, D'^'F+^, i,i = 2,3,..., (1.12) 

where D± = n_|_ • D and D± = n± ■ D. The fields in (jl.l2p carry twist £,£, and conformal spin 
1 -|- £/2, 1 + £/2, respectively, and their one-loop dispersion relations are given by (11.61 ll-7p . after 
plugging the appropriate value of the conformal spin s. For instance, the mass of the twist-£ excitation, 
equivalently £, is given to one loop by 

me = £ + V (V' (1 + em - V'(l)) + 0(/) . (1.13) 

The observation of |43j was further analyzed in [20j from which one can already extract an all-loop 
representation for the masses of the fields (|1.12p . Extending this analysis, we will construct the all- 
loop dispersion relations for the fields (I1.12p and argue on top of them that the excitations (I1.12p 



I am grateful to Juan Maldacena for pointing tliis reference out to me. 
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are bound states of gauge fields. One can already verify at the one- loop level, by plotting (jl.lSp . 
that rrii < ^rrigauge field = irn£=i, which is obeyed for a bound state. The relation between higher- 
twist excitations and gauge fields is nevertheless obscured at weak coupling because the i gauge-field 
constituents of the twist-^ excitation do not all live on the same rapidity sheet (idem for i by 0(2) 
symmetry). It becomes however manifest as the coupling gets larger and within a suitable kinematical 
range where the dispersion relations for bound states can be obtained by fusion of the dispersion 
relation for gauge field, in analogy with bound states of magnons over the BMN vacuum [44j though 
not in the same kinematical domain^ 

This completes the pattern of excitations considered in this paper. They are the momentum carriers 
for the large spin effective Bethe ansatz equations. There are also isotopic roots associated with the 
SO (6) symmetry of the GKP vacuum. Their embedding into the asymptotic Bethe ansatz equations 
is asymmetric because the BMN vacuum breaks spontaneously S0{6) down to SU{2) x SU{2). The 
root associated with the broken SU{2) C 50(6) is embedded in the form of a stack [33] which becomes 
isotopic at large spin and restores the symmetry. This mechanism is quite similar to what happens 
for spinous in the thermodynamic limit of the antiferromagnetic Heisenberg spin chain, where the 
restoration of the SU (2) symmetry is implemented by the formation of 2-strings of Bethe roots. 
For the stacks previously mentioned we verify that they do not carry energy nor momentum at any 
coupling. We show moreover that their equations are identical to the middle-node equations of the 
inhomogeneous SO{6) Heisenberg spin chain, with length fixed by the total number of scalars and 
with inhomogeneities given by the scalars' rapidities. 

Finally, let us comment on some underlying assumptions. Our analysis relies on the all-loop Bethe 
ansatz equations O [3] . It is well known that these equations do not capture correctly the interactions 
that wrap around the spin chain [50]. For a generic state, the size of these corrections is correlated 
to the spin-chain length and becomes more important as the length gets short. It seems therefore 
that our analysis is necessarily bounded to long operators, which is not welcome as we only assume 
that the spin is large. However, it is widely believed that wrapping issues are circumvented in the 
large spin limit. This is supported by computations of wrapping corrections of twist L anomalous 
dimensions in the s[(2) sector, using the Liischer formalism |5H I52j . These results indicate that 
wrapping corrections are delayed at large spin, 5 ^> 1, affecting contributions of order 0(log^ S/S"^). 
From our point of view, these contributions are effectively of a wrapping type as they are exponentially 
suppressed in the physical length R = 2 log 5. Hence we expect that the asymptotic Bethe ansatz 
equations account properly for the anomalous dimensions of large spin operators as long as we restrict 
ourselves to contributions suppressed by inverse powers of log S. It remains, nevertheless, conjectural 
that wrapping effects are suppressed by powers of the spin for all large spin operators. A more careful 
analysis would entail the use of the Y system of which is beyond the scope of this paper. 

The plan of the paper is the following. In Section [21 we present the embedding of the various 
excitations into the asymptotic Bethe ansatz equations [2j. We comment on subtleties in the sep- 
aration between isotopic and momentum-carrying roots at one loop, which are tied to the peculiar 
embedding of fermions. In Section O we perform the large spin analysis of the all-loop Bethe ansatz 
equations applying the technology of [lEl [3l |19] . In Section HI we derive exact representation for 
the dispersion relations of the various excitations. Then, we analyze these dispersion relations at 
weak coupling and derive explicit expressions for them up to three-loop accuracy. We discuss their 
analytical properties in the complex rapidity plane and derive fusion relations for bound states. In 
Section [SJ we perform the analysis at strong coupling and identify the various kinematical regimes. 
Section [6] contains concluding remarks and outlooks. Several Appendices contain the technical details 
of the analysis. 



^*The fusion is performed here by shifting the gauge-fields rapidities in the complex plane assuming u'^ < {2g)^, with 
u the center-of-mass rapidity (taken real), while it is done for > {2g)^ for bound states of magnons. 
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2 String hypothesis for GKP string 



We start recalling the asymptotic Bethe ansatz equations and then detail the embedding of the various 
excitations. We deal separately with the momentum-carrying roots and isotopic roots. We further 
discuss some subtleties in the embedding of fermions, that we relate to an enhancement of symmetry 
from SU{4:) to S'L(2|4) at weak coupling. 



2.1 Beisert-Staudacher equations 

The asymptotic Bethe ansatz equations exist in four distinct forms corresponding to different choices 
for the Dynkin diagram of the = 4 superconformal algebra p]. They are labelled by the gradings 
^1)^2 = il among which we select the one with rji = rj2 = —1 for our analysis. With the latter 
choice, the momentum-carrying roots are directly associated with the light-cone covariant derivatives 
[2], with =1)^1 in the bi-spinor notation of [10]. The equations in this grading are given 
by [2] 

l-Y[ ^2,fc - U2,j + i -Q U2,k - U3,j - 2 -Q ^2,fc - Uij - ^ 
^2 „, „, i K4, — 

I — TJ ~ ~ ^ TT ~ 

7=1 ^3,fe - U2,j + § jJi Xs^k - X^j 

9^/x, 

2 / 

''4,fc/ jjtk ■^4,fc •^4,i ^ ~ 9 /^4,A;-^4,j 



U2,k - U2,j - I .^^ U2,k - U3J + 2 j=l ^2,fe " Uij 



L K , — +1 9/4-- 

1 -r-r 4 k 4 1 hi / -^A '-•^'^ 



n"-4,fc ■^4,j y / -^4,^-^4 J 2/ \ 
^ —r. 2/ - + ^ Kfc,W4,i) 

-''/I I. •''4 4 J- U I X^ 

(2.1) 



^ n n n 4,^"'^^^ j-j- ^ 

j=i ^4,fc ~ ^3,i a^^^fc — x^j .^-^ 1 — g^ lx^-^x\^j 1 — g'^/x^^f.xjj 

-TT ^6,fc - U6,j + ^ -TT ^6,fc - ^5,j - 2 TT ^6,fc " '"TJ " 2 
„, „, i K4 T _ „2 



1 



U7,k - u(i,j - i -p4 1 - gVX7,kX, 



-TT "-/.fc "-b.j 2 TT _ 

M M7i.-Ufi , + i 1 



4,i 



jVfc ^7,fc - '^ej + § J=i 1 - g'^/xT^kxtj 
where x^ = x^{u) = x(n it |) and with x{u) the deformed spectral parameter 

x = \{u+^u^-{2gY), u = x+^-, (2.2) 

with asymptotics x ~ w at it ~ cx). The phase (t(u, v) entering the middle-node equations for the main 
(momentum-carrying) roots u^^^k is the BES/BHL dressing phase OES]. Its expression is recalled in 
Appendix |Al At one-loop, 5 — )• 0, we have both a{u, t;) — >■ and x{u) = u, and one recovers the 
Bethe ansatz equations for the psu(2,2|4) Heisenberg (super) spin chain [9]. 

Once a solution to Eqs. (|2.ip is known, the anomalous dimension of the associated operator is found 
as [MIE] 

5^ = 2g'Y:{^-^] , (2.3) 

j=l \^4,fc ^4,fe/ 
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where only the main roots (the spin-chain momentum-carrying roots) enter. More generally, the total 
spin-chain higher conserved charges are given by 



Qr+l = -i:(^-^] ■ (2.4) 

The global charges for a Bethe state with excitations numbers (i^i, . . . , K^) are given in [2]. The 
states we aim at considering are characterized by a large number of main roots ^ 1, but with 
Kj = 0(1) for j 7^ 4. All these states carry a large Lorentz (or total) spin 5 = ^ (si + S2) = K4+0{1), 
where (si, S2) are the Dynkin labels (twice the spins) for the su(2)^ C 5u(2, 2) symmetry algebra [2]. 
The charges that are kept of order 0(1) in the large spin limit include the twist, 

twist = Aq- S = L + Ki+ Kj . (2.5) 

where Aq is the canonical dimension. Another such charge is the transverse so(2) spin = ^(s2 — si) = 
^ {Ki + — Kj — K^). Of course, all the su(4) quantum numbers are of order 0(1) in this limit, 
and the GKP vacuum is a singlet under this symmetry. The most important observable is the GKP 
energy given by 

E = A- S = L + Ki + K7 + 6A, (2.6) 

where 5 A = 0{g^) is the anomalous dimension. Comparing (j2.5p and (j2.6p , we see that the twist 
coincides with the energy of a state at tree level in the gauge theory. However, starting at one loop, 
the total energy is no longer a finite quantity in the large spin limit. The energy measured above the 
vacuum is however finite and reads 



E - -Evacuum = twist - 2 + 5A- (5Atwist-two , (2.7) 

where 

^vacuum = 2 + (JAt^^-two = ^2(5) (log S - ^(1)) + 2 + ^2(5) + 0(log S/S) , (2.8) 

with ^2(5) = 2rcusp(5')) V'(l) = ~7e is the Euler-Mascheroni constant, and finally -62(5) = 0{g'^) is 
the so-called (twist-two) virtual scaling function |53j . The expression (|2.8|) is slightly cumbersome, 
but it will not be difficult to make use of it. This is because, in the large spin limit, 5^1, the 
Bethe ansatz equations (|2.1|) can be turned into a linear integral equation for the distribution density 
of roots [351 El [El Is]. Thanks to the linearity of this equation, the separation of the various 
contributions to the energy in ()2.7p becomes straightforward. 



2.2 Momentum carriers 

The embedding of the momentum carriers, which were presented in Introduction, can be obtained 
from the literature [17^ [T9l [Ml [20] on the large spin limit of the equations (I2.ipl^^l The scalars, for 
instance, are well-known to be associated with the dual solutions of the main roots equations [17| and 
they appear as holes in the distribution of magnons mode numbers [19] . The total number of holes 
in a state counts the total number of scalars on top of the GKP vacuum!^ One can easily introduce 
fermions and ^+ by means of the roots ^3 and us, respectively. Indeed, adding a root M3, for 
instance, performs the change of flavor D^Z — t- This is not yet the end of the story for the 

embedding of fermions and there are subtleties that we discuss below. The embedding of the gauge 
fields and of the higher- twist excitations (I1.12p was analyzed in [20] . It was found that they come 
out as stacks of strings, that are complex solutions to (j2.ip composed of roots of various flavors [3S]. 
For instance, the gauge field with 50 (2) spin 1, is embedded in the form of a 2-string of roots 

Us centered around a real root U2. Similarly for F^± with the replacement U3 — t- U5 and U2 — ?■ uq. 

^^Except for the fermions, whose embedding turns out to be subtle, as explained in a separate subsection. 
^^Note that we do not take into account the two large holes [171 119j which cannot change position in the lattice of 
mode numbers. They form part of the background. 
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Of course the strings that compose the stacks are not 'exact strings' at finite spin. But in the large 
spin limit, they are, as observed in [20] and as we would expect from analogy with thermodynamical 
limit of compact spin chain. The estimate for the departure from the exact string is expected to 
be exponentially suppressed with the (physical) length, i.e. associated to ~ 1/5 corrections. Now 
the gauge fields are just the first representants of two infinite series of higher-twist excitations |43] . 
associated in (|1.12p to partons of the type D^^^K^x. From the analysis of [2D] we deduce that the 
twist-£ excitations is embedded as a stack of a {i+ l)-string of roots U3, a ^-string of roots U2, and a 
{i — l)-string of roots Ui. Similarly for after the wing exchange (ni,n2,?i3) — ^ (nyjUgj^s). 

As already mentioned in Introduction, we interpret the twist-£ excitation ~ as a bound 

state of £ gauge fields. Their form is however a bit puzzling and not immediately evocative of their 
interpretation. For a bound state of i gauge fields, one would have probably expected something like 
{F+±Y. There seems however to be no room in the ABA equations for accommodating the latter 
parton as an excitation on the large spin background. Indeed, looking for such an excitation, we 
would turn to the (s = 3/2) s[(2) spin chain Ti D^F^^ and ask for a solution describing the bound 
state ~ {F^±Y on the twist-two background - here realized by ~ Tr D^F^j_. The sought solution 
should then be given by a distribution of £ hole^^ that is expected to form a string in the complex 
rapidity plane, in analogy with bound states of Bethe rootscj. But no such complex solutions exist 
for the non-compact spin chain [35j. The operator D^j^^F^^, which is not contained in the latter 
subsector but in a larger 5u(2, 1) subsector that supports string solutions, appears as a substitute for 
the 'bound state' {F^±Y in the spin-chain description. 

At the same time, how should we understand the claim that the operator D^j^^F^± is a stable 
excitation over the large spin background. Naively, if we form a spin chain out of such operators, one 
would certainly expect that the covariant derivatives D± would hop from one site F^± to another, 
suggesting that the numbers of bound states are not separately conserved. This is not very welcome 
as, if the large spin limit is associated to an integrable model of sort, the number of excitations with 
a given mass (and quantum numbers) should be well defined. The mixing alluded to before certainly 
occurs, but the question is whether it persists at large spin and for the asymptotic part of the wave 
function, when the excitations are separated from one another by a large number of background 
derivatives D^. Without the explicit form of the large spin wave function for such states, it is 
difficult to rigorously address this issue. We will give nevertheless a somewhat dynamical argument 
in favor of this interpretation of the twist-^ excitation as a bound state of gauge fields, later on. It 
will be based on the fusion identities of the type 

e-i 

Et{u) = E^fiu - i(l - l)/2 + ij) , (2.9) 
i=o 

and similarly for the momentum, which are natural for bound states. These identities will be observed 
to hold only when < (^gf' and thus will not reveal themselves at the level of the weak coupling 
expansion, which always requires > (2^)^. The obstructions met in extending these relations at 
weak coupling are tied to the presence of branch points in the complex u-plane, that makes the shifts 
along the imaginary n-axis in (j2.9p inequivalent if performed inside or outside the interval < g^. 
This feature might explain the complications mentioned above in the realization of the bound states 
of gauge fields in the spin-chain picture, since the gauge-field rapidities in a given bound state do not 
all live on the same rapidity sheet for u'^ > {2gY . 

2.3 Isotopic roots: one loop versus all loops 

So far we have been considering the roots ui,U2, uq,ui, as elements in the construction of stacks. Now 
we would like to examine them separately. We naively expect that large-spin solutions to the Bethe 

^'^Note that the holes introduced here are distinct, though very similar, to the ones associated to the embedding of 
the scalar excitations discussed before. 

'^^To form a bound state we need to find a way to store less energy in the state than allowed by a configuration of £ 
real holes - hence the need to put some holes in the complex plane. 
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ansatz equations ()2.ip with solitary (real) roots ui, U2, ug, uy, do exist, and so we need to interpret 
them. We start with the roots U2,uq, whose meaning is quite clear. 

The roots U2,uq, play exactly the same role for the GKP operators as they do for the BMN ones. 
Namely, they are associated with the SU{2) x SU{2) C S'C/(4) (residual) symmetry group that is 
common to both backgrounds. They are isotopic roots that do not carry energy nor momentum and 
implement the latter two SU (2) symmetries. But the residual symmetry of the GKP vacuum is SU{4) 
not SU{2) X SU{2). We expect therefore a third type of isotopic root. The missing root is embedded 
in the form of a stack of roots M3, tt4, 115. More precisely, the stack is made out of one 2-string of main 
roots, U4^± = Mb =t i/2, with one root U3 and one root M5 on top of each other, U3 = = Ub- It is 
not obvious that this stack does not carry energy nor momentum. This is true however to all loops, 
as shown in Appendix O What happens is that the bare energy, carried by the 2-string, is exactly 
compensated by the induced energy, which measures the shift in the energy stored in the large-spin 
continuum of real main roots due to the presence of the stack. Moreover, the large-spin equations for 
these stacks take exactly the same form to all loops 

i K\, . i ^6 i 

-r-r U^,k - Uh,j + 2 ^ TT ^b,fc - U]^,j + i TT Uh,k - U2,j - 2 TT ^b,fc - Uqj - ^ 

jJl Uh,k - Uh,j - f f^_^j Ub,k - UbJ - i .^-^ U],^k - U2,j + I Uh,k - Uqj + § 

as shown in Appendix O Here is a hole (scalar) rapidity, and K2,Kq, are the numbers of solitary 
roots U2,uq. There is no explicit dependence on the spin in these equation, and so there is no momen- 
tum for a root u^. We recognize in (|2.1U|) the central-node equations for the inhomogeneous SO (6) 
Heisenberg spin chain, with length equal to the number of scalars and with the scalars' rapidities en- 
tering as inhomogeneities0 We conclude that the stack Ub effectively restores the SU{2) symmetry, 
initially broken by the BMN vacuum^ In essence, the mechanism that restores the SU{2) symmetry 
here is identical to the one that operates in the thermodynamical limit of the anti-ferromagnetic 
SU{2) spin chain. The only difference is the necessity here to have fermionic roots to support the 
2-strings. 

The SU{4:) roots U2,uq, and u^, are isotopic to all loops. But they do not form the complete 
set of isotopic roots at one loop, however. Indeed, at one loop, the fermionic roots ui,uj, are also 
isotopic. They do not carry momentum and their energy is canonical and equal to 1, with no one-loop 
correction. The reason being that they do not couple to the middle-node equations, at this order, and 
so cannot induce contributions to the energy. These extra isotopic roots enhance the symmetry from 
SU{4:) to S'L(2|4). Under this symmetry, each momentum carriers give rise to an infinite dimensional 
multiplet, with the 'infinite direction' coming from the possibility to extract an arbitrary number 
of covariant derivatives D_. At higher loops, however, roots ui,U7, couple to the main roots, see 
Eqs. (j2.1|) . and thus acquire energy and momentum. It raises the question of their role in the spectum 
of excitations. As we shall see they are in fact fermions, which are trapped in a small momentum 
domain p = 0{g^) at weak coupling. 

2.4 Fermionic subtleties 

The fact that the roots ui,uj, are fermions is supported by the observation that these roots can actu- 
ally act on a covariant derivative at higher loops, due to the supersymmetry transformations [10] 



5iD+^g-^+, 6rD+^g^+. (2.11) 



^^Note that the scalars' rapidities are the only inhomogeneities in the central-node equations (|2.1U[) . in agreement 
with the fact that the scalars are the only excitations in the 6 of SO{6). Fermions, which are in the 4 or 4, couple to 
the equations for roots U2 and ue, and gauge fields and bound states, which are neutral, do not couple to any of the 
SO{6) isotopic roots. 

^°Note that at one loop, we could as well work over the pseudo- vacuum Tr_F^j^, associated to the choice of the 
distinguished Dynkin diagram (the Beast of 0). The SC/(4) symmetry would be manifest from the very beginning. 
It is not known however how to do it at higher loops. Similarly the choice of the pseudo-vacuum Tr would keep a 
SU{3) C SU{4:) symmetry subgroup manifest. 
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Accordingly, the one-loop 5L(2|4) descendants, obtained by acting on the momentum carriers with 
the roots ui,u-/, presumably decay into fermions at higher loops. 

At the kinematical level, it happens that the physical rapidity is not the Bethe rapidity but its 
Zhukowsky map x. The x-plane requires to glue together two u-planes, and one finds two expressions 
for energy and momentum of a fermion with rapidity u, one coming from ui and the other one from 
Us, idem with uj and U5. For illustration, we find 



£:,f(n) = l + %^ + 0(/) 



Psf{u) 



u 



for a small fermion with < and 



Euiu) = 1 + 2^2 (^(1 + iu) + - iu) - 2V'(1)) + 0{g^) , pif{u) = 2u + 0{g^) 



(2.12) 



(2.13) 



for a large fermion with > 9^- These two cases are realized at weak coupling by the roots ui and 
U3, respectively, with x ^ g"^ /ui in the first case and x ~ ua in the second one. The two are however 
connected with each other, see Fig. [2J For instance, in both cases the dispersion relation reads 



E{p) = 1 + 2/ (^-(1 + ip/2) + V(l - ipl2) - 2^(1)) + 0(/) , 



(2.14) 



but while p = 0(1) at u fixed and 5 ~ for a large fermion, it is small p = 0{g^) for a small fermion 
under the same assumption. In particular, a fermion can be found at rest if it lies at x = 0. Then its 
mass is exactly one 

ruf = Ei{p = 0) = E{{x = 0) , (2.15) 

due to the special property of the point x = in the all loop asymptotic equations [5]. At one loop 
this point can be identified with u = 0, but at one- loop only. 



fermion at rest 



small fermion brancii 
X < g 





" = 91 






= 


\u = 2g 


u > 2g 


u - 



X > g 
large fermion branch 



quasiclassical limit 



Figure 2: The two fermion branches. 

For some dynamical reason, tied to the coupling with the SU{2) x SU{2) roots [2j, some fermions 
get trapped at weak coupling in the small momentum domain and then behave essentially like if they 
were isotopic. The separation between roots 1*3,^5, and ui,U7, seems therefore to be a dynamical 
necessity at weak coupling. However, for our purpose of constructing dispersion relation at finite 
value of the coupling, it is quite artificial and it is better to restore the unity between fermions. This 
can be done by using the dynamic transformations of [2J . The price to pay is that we need to change 
the length L in (|2.1|) as 

L + K[ + Kt, (2.16) 
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where K[ and K'-j are the numbers of (sohtary) roots Ui,Ut, transformed into roots u^jU^. The 
quantity above counts the number of length-one partons in a given states and it is equal to 2 + M 
where M is the total number of momentum carriers and 2 comes from the two vacuum-building 
partons. 

These remarks about the fermions complete our discussion of excitations over the GKP vacuum. We 
are still lacking a proof of the completeness of the basis presented before. However, it fits remarkably 
well with the analysis of [56], which deals with the string hypothesis for generic super spin chain, and 
with the independent analysis of the spectrum around the GKP string performed in [21] . 



3 All-loop integral equation 

Assuming the pattern of excitations presented in Introduction, and the embedding discussed in Sec- 
tion [21 the next step is to derive the dispersion relations. This is done by analyzing the solution 
to the middle-node equations in (12. ip for a large number K4 of main roots and for a generic set of 
excitations. This task is greatly simplified by the large spin limit [351 [T7] in which the latter equations 
can be turned into a linear integral equation for the distribution density of main roots [181 US] • Our 
analysis, in this section, follows the methodology of [181 |3l [19] to constructing this equation to all 
loops. We perform a few generalizations, as we need to treat both the symmetric and antisymmetric 
part of the distribution density, and simultaneously include all the excitations encountered previously. 
Incidentally, we observe that the symmetric and antisymmetric part totally decouple from one an- 
other, leading to integral equations of the type already considered in the literature [18 1 13[ [T9 | [53 [ 120]. 
Applying to them the tools developed in [211 [Ml [Ml [Ml [33 [M] for dealing with these equations, 
we derive dispersion relations for the various excitations. Our construction does not provide explicit 
analytical expressions for the all-loop dispersion relations. But it gives compact expressions for them, 
parameterized in terms of the solution to the BES equation, which controls the vacuum distribution 
of roots to leading order at large spin. A sample of the literature which is relevant for the analysis 
presented below includes [351 [T71 [THl [TUl [2D] and the reviews [57] . 

Some preliminary remarks are in order. In the following, we assume that the dynamic transfor- 
mations [2] have been performed on the solitary roots ui,uj, to unit the fermionic roots. These 
transformations lead to the replacement of the spin-chain length by 

00 00 

L + K[ + Kl^ = 2 + M , M = Ki, + Ki + Kj + Y,Ki + J2Ke> (3-1) 

e=i i=i 

where K[,K!f, are the total number of solitary roots ui,ui, respectively, M is the total number of 

momentum carriers, and is the total number of scalars (* = h), fermions (* = f), There 

is actually no need to separate between excitations carrying positive or negative 50 (2) spin. Since 
the main roots do not distinguish between them, they are degenerate, and it is enough to consider 
★ = h,i,i. Finally, the particular case of the SU{2) isotopic root will not be treated below. Its 
analysis is deferred to Appendix [Cl where it is shown that it does not carry energy/momentum, nor 
any of the spin-chain conserved charges, to all loopsl^ 

3.1 Middle-node equations 

In the presence of a generic set of excitations, the middle-node equations in (|2.1|) can be written 

f \ "^^^'^ ^4 M 

1= ^^ J]^ 5|4(n4,fc,M4j)cr^(M4,fc,M4j) 54^(U4,A;,'U^) , (3.2) 

^^There is a slight contribution to the total spin-chain momentum e'^ in the form of the factor (— where K^o is 
the total number of 5(7(2) isotopic roots. 
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where all the ^4 main roots are realo For later convenience, we have factorized the scattering phase 
for the self-interaction of main roots into its undressed component, denoted Sl^{ui^k,U4j), and the 
dressing phase, 0"^(u4^fc, U4j). The symbol -k denotes a generic excitation, which can be a scalar (hole) 
* = h, a fermion ★ = f , a gauge field * = £ = 1 or a bound state thereof -k = i > 1. The total number 
of excitations within the state is fixed equal to M and the phase 54* (1*4, u^) stands for the scattering 
phase between a main root U4 and an excitation -k with rapidity u^. The product over the excitations 
is restricted to ★ 7^ h because there is no scattering phase for scalars at this time. They are indeed dual 
to the main roots [17j and appear as holes in their distribution of mode numbers p!9]. A scattering 
phase for them will be induced later, during the process of taking the large spin (continuum) limit. 
For any of the other excitations, the scattering phase with a main root is obtained directly from the 
middle-node equations in (j2.ip . taking into account the embedding of each excitation into the Bethe 
ansatz equations. For gauge field {i = 1) and bound states {i > 1), which appear as stacks of various 
strings, they are obtained by fusion of the relevant elementary scattering phases (see Appendix IB] for 
details) . 

Now we would like to analyze the equations (|3.2p at large spin, S = K/^ + 0(1) » 1. Following the 
lines of [181 [T9] , we cast them into the form 

2Tmk = Z{uk) , (3.3) 

which follows from taking the logarithm of (|3.2|) . after multiplying both sides by the phase (— l)^4+i^h-i 
leading to a more convenient choice for the branch of the logarithm. Then the (fermionic) mode num- 
bers Uk, entering (j3.3p . belong to a finite lattice (3 of integers or half-integers - depending on the 
parity of + — 1 - and the so-called counting function Z{u) is defined by 

Z{u) = Z{0) + 27r r dv p{v) = -iy log (-5^4 (u, U4,,)) + . . . , (3.4) 

where the dots includes contributions from dressing, that we will restore later. Meanwhile we 
have introduced the density distribution p{u) as the derivative, with respect to the rapidity n, of the 
counting function. The lattice of mode numbers & is almost filled up by the magnons up to a few 
holes [19j. These holes are associated to dual solutions [17J of the equations ()3.3p and they can be put 
in correspondence with a set of rapidities by using (j3.3p . Namely, a hole occupying the mode number 
nh in S is given the rapidity Uh = Z~^{2TTn\^). We are assuming here that the counting function 
Z{u) can be inverted, at least for real rapidity n, and thus that the mapping between rapidities and 
mode numbers is one-to-one. This is definitely the case in the s[(2) sector, which is obtained by 
imposing the restriction that the only momentum carriers are the holes (scalars), M = K\i with i^h 
the total number of holes. This is not obviously the case here, but, as we shall see below, this is a 
quite reasonable assumption at large spin. We will therefore push further the analogy with the s[(2) 
sector and now take the large spin limit. 

The main simplification at large spin limit is that the equations (13. 3p can be turned into an integral 
equation for the distribution density p{u). The first step in the derivation of this equation is to 
approximate the sum over the main roots in (13. 4p by an integral weighted by the density. To this 
end, we extend the latter sum to the full lattice of mode numbers, subtracting simultaneously the 
holes' contribution, 

- i X log (-544 (u, UAj)) = log (-5*44 (u, u{n))) + i ^ log {-Sl^{u, Uhj)) , (3.5) 

i=i nee j=i 

where u{n) = Z~^(27rn) is the rapidity, either of a magnon or a hole, corresponding to the mode 
number n in (3. Now in the large spin limit, the distribution of roots u{n) is dense [35^ 117] . the 



^^In our analysis, complex roots 114 appear in the form of 2-strings, which are part of the 5*17(2) isotopic roots, 
considered separately in Appendix [Cl 
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typical distance between the roots of order 0(5") being ~ l/log5. It allows us to apply the Euler- 
MacLaurin summation formula to evaluate the first sum in the right-hand side of (|3.5p . giving 




where dots stand for 1/5-suppressed corrections at large spin. Note that we have assumed that 
the distribution density p{u), as defined in (I3.4p is positive over the full domain of integration. This 
property, which is related to the invertibility of the counting function alluded before, is not guaranteed, 
but seems quite realistic, in the large spin limit. The reason being that the distribution density p{u) 
receives a dominant contribution from the vacuum (twist- two) distribution of roots. The latter scales 
like [351 [El [IS] 

/5vacuum(n) = - log 5 + 0(logO S) , (3.7) 
TT 

at large spin, for a rapidity u = 0{S^), which is the regime of interest here. When holes, or other 
excitations, are introduced with rapidities of order 0{S^), they generate corrections to the distribution 
density p{u) which are of order 0(log° S), and thus subleading as compared to ()3.7p . The density, as 
defined in (j3.4p . should therefore be positive, at least for u = 0{S^). 

It remains to specify the domain of integration in ()3.6p . For our purpose, it can be considered 
to be the full real line, but it requires some comments. Indeed, it is known [3S] that the twist- 
two distribution is supported on the symmetric interval < (5/2)^, at large spin. It does not 
immediately implies, however, that we can extend the domain of integration in (|3.6p to ±00, even if 
we are interested in constructing the density p{u) for u = 0{S^), which is our ambition here. The 
reason is that the vacuum density scales like 

Pv.cunm{u) = - log (S/u) + 0{l/u) , (3.8) 
TT 

at large rapidity. It is so to all loops [18j, as long as S" ^ 1 and S ^ g, which are both always 
implicitely assumed here. Then one needs to proceed with care concerning the domain of integration. 
But, fortunately, the situation is different for corrections to the density sourced by the excitations. 
Indeed, the latters are suppressed at large rapidity, if the excitations carry rapidities of order 0{S^). 
It follows that in our case, we can effectively extend the domain of integration to the full real line. 
We note finally that this implies that the large rapidity regime of the density p{u) is controlled by 
the vacuum contribution and scales like (j3.8p . One can think therefore of the behavior (j3.8p as a 
boundary condition for our problem. One can change this boundary condition by sending a few 
excitations (holes, ...) to infinity, or more precisely by giving them rapidities of order 0{S). The 
boundary conditon (|3.8p then gets replaced by [3S] 

(2 + n) 

p{u) = ^— ^ log (S/u) + Cn + 0{u) , (3.9) 

TT 

where n is the number of large rapidity excitations and where Cn is some constant. This corresponds 
to working around higher-excited trajectories in the spectrum of large spin scaling dimensions [16^117). 
A recent literature on this subject with analysis performed at both weak and strong coupling, and 
with relation to spiky strings [58| , includes [361 EHl EQ] . 

The next step in the derivation of the integral equation for the density p(u) is to differentiate with 
respect to u in ()3.4p . using ()3.6p to evaluate the right-hand side of the last equalities in ()3.4p . This 
step is performed in the next subsection after restoring contributions associated to dressing phase, 
.... in (f^a . 

3.2 Large spin integral equation 

The complete counting function, deduced from the equations (j3.2p . can be decomposed as 

Z{u) = Z°{u) + {2 + M)Zicngth{u) + ^dressmg(w) + ^excitations («) , (3.10) 
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where the first term, Z°{u), coincides with the one previously considered in ()3.4p . It is the term that 
generates the main kernel of the (large spin) integral equation |18[ [3] . Its expression follows from the 
relation [TBI 



zlogi-Sl,iu,v)) = -^log - 2zlog ( ^^^/^] , (3.11) 



1 — iu + ivJ \l — g'^/x 

where x{u) = (n + y^u'^ — {2g)'^)/2,y{v) = {v + \/v^ — {2g)'^)/2 and = x{u ± i/2) (similarly for 
y^). The second term in (|3.1Up takes care of the contribution proportional to the spin-chain length, 

(X ( w) \ 
~^+(^ J • ^^'^^^ 

Concerning the contribution originating from the dressing phase, it is better not to expand it for the 
time being. So we simply write 

-^dressing (-u) = 2 ^ 9{u, U^j) , (3.13) 

i=i 

where 6{u, v) is given in terms of the BES/BHL dressing phase a{u, v) as 6{u, v) = —i \oga{u, v) . 

Finally, we have the contributions from the various excitations 

M 

-^excitations (li) = X! {-Sii,{u, U^)) . (3.14) 

The integral equation for the density is linear and, as we have said before, we are only interested 
in the contributions sourced by the excitations (holes, ...). We therefore decompose p{u) into two 
part@ 

p{u) = /9vacuum(w) " Cf{u) - ^{u) , (3.15) 

where the first term stands for the twist- two (vacuum) distribution of roots. It is controlled to leading 
order at large S by the BES equation [THl |3] and to next-to-leading order by a similar equation [53] . 
The other terms in (|3.15|) account for the deformation of the vacuum distribution induced by the 
excitations. The two densities (j{u) and ^{u) are univocally distinguished by their parity in u. The 
former is even, cr(— u) = (t{u), while the latter is odd, (j{—u) = —a{u). This separation anticipates on 
the fact that the equation we will have to consider is not of difference form. Now, implementing the 
decomposition (|3.15|) consistently requires to remove the terms in the counting function (|3.1U|) that are 
captured by the vacuum distribution density. This is done by performing the replacement 2 + M ^ M 
in (j3.1Up and by subtracting the vacuum contribution to the dressing-phase correction .^dressing (^^)- 
The latter quantity admits an expansion over the spin-chain higher conserved charges [2 [55l [3]. 
Substracting the vacuum contribution amounts to evaluating these charges above the vacuum, 

Q, = g;— + (3.16) 

where the charges Qr are the total charges induced by the excitations. In the large spin limit, they 
are found as 

qr = - du (o-(u) + a{u))qr{u) - ^ qriu\,j) , (3.17) 

i=i 

where qr{u) is the charge carried by a magnon with rapidities u. 



^The minus signs are conventional |18l 119) . 
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The equation for the distribution densities o"(u), a'{u), fohows from differentiating (13.40 with respect 
to n, applying (|3.1U|) and (|3.6p . Keeping only the relevant contributions for the analysis of the 
excitations (see paragraph above), we get 

= 2TTa{u) + 2Tra{u) + duZ{u), (3.18) 

where Z{u) is given by (|3.10p with the modifications explained before. The contribution to ()3.18p 
from Z°{u) splits into two pieces, after using ()3.6p . 

duZ°iu)\^^^^^ „^ = i dvdu log {-Sl,{u, v)) (aiv) + a{v)) + i ^ a„ log i-Sl,{u, Uh,,)) . (3.19) 

i=i 

The first term generates the so-called main kernel [3l [18] for the integral equation while the second 
one is the source term associated with the holes rapidities. Now, remarkably enough, the parity odd 
and even densities decouple from each other at the level of the kernel. Namely, we have 

i dvdu log {-Sl^{u, v)) {a{v) + a{v)) \^^^^ ^ = /C° * a{u) , 

_ (3.20) 
ij dvdu log {-S2i{u, v)) {a{v) + a{v)) \^^^ ^ = * a{u) , 

where explicit expressions for the two kernels /C°,/C°, will be given later. It is not clear to us whether 
this property should have been expected on general grounds. One could still expect a coupling at 
the level of the dressing-phase correction 9n-^dressmg('w), but here again there is none. The even part 
depends on the even charges and the odd part on the odd ones. So the equation ()3.18p splits into 
two integral equations that read 

= 27r o-(n) + /C° * a(u) + I(u) , 

(3.21) 

= 2TTa{u) + IC° *a{u)+I{u), 

where the inhomogeneous contributions I{u),Z{u), are respectively even and odd in u. Note that our 
terminology here is slightly abusive since the inhomogeneous terms contain the contributions from the 
dressing phase, that will correct the kernel of the integral equation later on. It is convenient never- 
theless to keep them enclosed in I(u),I{u), together with the genuine inhomogeneous contributions, 
at this moment. 

We are now in position to detail each term in the equation (j3.2ip . We start with the kernels. From 
their definitions (I3.20p . and taking into account the expression ()3.1ip and formulae in Appendix [Aj 
we get 



IC°-ka{u)= —2 J dt cos {ut)e * j dv a{v) cos (vt) 

+ 2{2gf / dt cos (ut) e-*/^ / ds tK{2gt, 2gs) e''^^ / dv a{v) cos (vs) , 

i° . ^° ^ (3.22) 

iC-ka{u)= -2 J dt sin {ut)e ^ j dv a{v)sm{vt) 

+ 2{2gf / (it sin (nt) e-*/2 / dstK{2gt,2gs)e-^/'^ / dv a{v) sin (vs) , 
Jo Jo J 

where K{t, s) is the symmetric, K{t, s) = K{s, t), kernel density of [TS1[3]. It can be decomposed into 
even and odd parts in t as 

K{t,s) = K+{t,s) + K^{t,s), (3.23) 
and the two components K-i-(t,s) can be found to admit the representations [18^ [3] 

K+{t, s) = ^Y.^2n- l)J2n-l{t)J2n-l{s) , 

2 (3.24) 

K-{t,s) = -J2{2n)J2n{t)J2n{s), 

n>l 
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where J2n-i{t), J2n{t), stand for the parity odd, even, Bessel's functions, respectively. 
The inhomogeneous terms in (|3.21|) can be written as 



M 



^(^) = '^I*{u,U^) +Xdrcssing('u) , 
* 
M 

i{u) = '^i^{u,U^,) +2d: 



M (3.25) 



where the sums run over all the M excitations above the vacuum, including the holes. For any 
excitation with rapidity n^, we have 



X (u) \ 

x+{u) ) ' (3.26) 



I odd in u ' 

where for a hole, * = h, the scattering phase has been induced in Eq. ()3.19p . 

- idu log (-54h(u, Uh)) = +idu log {-Sl^{u, Uh)) . (3.27) 

Note that we have distributed the contribution coming from duZ\ength{u) among the M excitations. To 
proceed further with the inhomogeneous terms (|3.26|) we need explicit expressions for the scattering 
phases between main root and excitations. Unfortunately, we did not find a closed formula that 
would encompass all these contributions, and each type of excitation has to be treated more or less 
separately. However, they have in common that they can all be written as 



l-CO 

X^(u, Ui,) = —2 dt cos {ut) e~*/^ 7r^(t, «*) , 
Jo 

Xi,{u,Ui,) = —2 / sin (tit) e~ / ^^(t,^^). 



(3.28) 



The expressions for the functions 7r^(t, u^), ■7f^(t, u^), are given in Appendix |Bl The contribution 
coming from the dressing phase can be written as (see Appendix IB]) 



Xd™g(u) = 2(2g)2 dt cos (ut) e-*/2 H dstK^{2gt,2gs f f^^f^ 
Jo Jo e*— i 

Xdressing(t^) = 2(2^)2 H dt siu (ut) c'^/^ dstK+{2gt,2gs)^±p^ 
Jo Jo — i 

where we have introduced generating functions for the higher conserved charges as 

q.{t) = 2 5](2n - 1)(-1) V"~''?2n J2n-l(t) , 

n>l 

Q+(t) = 2 5](2n)(-l)V^2n+l J2n(t) • 



(3.29) 



(3.30) 



n>l 



Note that the subscripts '— /+' refer to the parity in t, and not to the parity of the charges that are 
generated. We recall that the charges qr{r > 2), entering (|3.3U|) , are evaluated above the vacuum. 
We can now verify the claim concerning the coupling between cr{u) and d'{u) at the level of the 
dressing-phase corrections. Indeed, the charges q2n,Q2n+i, only depend on a{u),a{u), respectively. 

Looking at the various contributions to the integral equations (I3.2ip suggests to switch to a Fourier 
representation. Following [18\ [3]. we introduce it as 

n{t) = e~*/2 / du a{u) cos (ut) , 0(t) = e'^/^ / du a{u) sin (ut) , (3.31) 
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M 



(3.34) 



for t > 0, and reciprocally 

1 1 r°° 

a{u) = - dt cos {ut) e*/2 9.{t) , a{u) = - dt sin {ut) e*/^ 0(t) . (3.32) 
vr Jo TT Jo 

Note that the functions J7(t),r2(f), do not have a definite parity in t. Taking the representations 
above into account, the integral equations ()3.2ip can be brought into the form 

0=/ dt COS {ut)f{t) , 0=/ dt sm{ut)f{t) , (3.33) 
Jo Jo 

with some functions f{t),f{t), that collect all the terms given before. The equations (j3.33|) are 
valid for any real value of u by construction. Under some minimal assumptions of regularity for the 
functions f{t), f{t), the unique solutions to (|3.33|) are f{t) = f{t) = for any t > 0. Expanding now 
the expressions for f{t) and f{t), we get the equations for il{t),Q{t), as 

POO 

{e^ -l)Q{t) + {2gf / dstK{2gt,2gs)n{s) 
Jo 

+ {2gf r dstK^{2gt,2gs)^^Pp- =J2'^At, 
Jo -1 

roo 

(e* -l)$7(t) + (2g)2 / dstK{2gt,2gs)n{s) 
Jo 

+ {2gf r dstK+i2gt,2gs)^±^ = f2ir.it, 
Jo e i 

To make the picture complete, we rewrite the generating functions for the conserved charges (|3.3U|) 
in terms of 0,{t), Q{t) ()3.3ip . Starting from ()3.17p and using formulae in Appendix[Al we get 

roc 

q-{2gt) = 2{2gf / dstK+{2gt,2gsMs) 
Jo 

+ 2(2^)2 J2 / ds tK+{2gt, 2gs) cos (uh^s) e'^/^ , 

.00 (3.35) 

q+{2gt) = 2{2gf / dstK^{2gt,2gs)n{s) 
Jo 

+ 2{2gf y / ds tK^ {2gt, 2gs) sin (uh js) e''/"^ . 

The equations (j3.34p are of type previously encountered in the literature [181 El [13 [53] . In the 
following subsection, we apply the methodology that has been developed to deal more conveniently 
with them in [6ll |62l [Ml [Ml [Ml [65] . 

3.3 Simplifying the equations 

The equations (|3.34p being linear both in the unknows 0,{t), Cl{t), and in the sources ^) 
can treat each contribution separately. Thus we look for a solution in the form 

M M 

n{t) = J2 > ^(t) = E ^*(^) ' (3-36) 

where r2^(i), r2^(t), are sourced by a single excitation carrying a rapidity n^. In parallel, we have a 
similar decomposition for the conserved charges, or equivalently their generating functions (|3.3Up . 

q. {2gt) = E 9- (29*) , 9+ (2^0 = E (2^*) , (3.37) 
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where in particular case of the holes, * = h, we need not to forget the extra contribution in the 
right-hand side of (|3.35|) . 

Now, to turn the equations for into a form which is more appropiate for our purposes, 

we introduce auxiliary functions 7j_(t), 7^(t), with definite parity in t, '~f±{t) = ziz'j^{—t) and similarly 
for 7±(i). This is done schematically as 

(e* -l)n^{t) = one-loop value + ■y'L{2gt) + 7|(2c/t) , 

, _ (3.38) 
(e* -l)n^{t) = one-loop value + 7* (2^*) + 7|(2c/t) . 

This is sketchy because, actually, it is convenient to add a bit of higher-loop corrections into the 
'one-loop values' above. The precise decompositions can be found in Appendix [Bl for each type of 
excitations (|B.5l IB. 151 IB. 241 IB.34|) . Plugging the ansatze (|3.38p into the integral equations ()3.34p , one 
finds that the functions 7^(t), 7j_(t), are naturally given in an expansion over the Bessel's functions, 

7:^(t) = 2^(2n-l)72Vi^2n-i(t), 

"-^ (3 39) 

7;(i) = 2 5](2n)72*„J2„_i(t), 

n>l 

and similarly for ^'±{t). The expansion coefficients 7^,7^^, are functions of the coupling constant g, 
and they implicitely depend on the rapidity n^. Note also that there is no contribution ~ Jo{t) for 
7+(i),7+(t). It implies that 7+(t) ~ 7+(*) ~ at small t. For 7*(i),7*(t), we have the smah t 
asymptotics 7* (t) = 7i((?)t + 0{t^), and similarly for 7i(t). 

Now, adopting the notations (|3.38l I3.39|) , the integral equations (|3.34p can be casted into the form 
of infinite systems of equations for 7^,7^, 



00 



Jo t e i , . 

7n + / —Jn{2gt) -— = K^(U^) , 

Jo t —1 

with n > 1. The two series of coefficients (u^), (u^), are given in Appendix [Bl for each type 
of excitations (|B.6l IB. 161 IB. 171 IB. 251 IB. 261 IB.35p . In this form, the equations are amenable to a 
numerical analysis along the lines of [61j. For our purposes, it is the relationship of (I3.40p with the 
anologue system of equations for the vacuum solution that is important. As we have said before, the 
vacuum distribution density is controlled, to leading order at large S, by the solution to the BES 
equation [3| . This equation is known to support the form p.40p , see [STl E21 ES] , with 

K^ = 2g6n,u ^n=0. (3.41) 

The two functions 7±(t) that solve (|3.4Up with the source terms (|3.4ip . are fundamental within our 
analysis 1^ Indeed, as we shall see, they control the dispersion relations of all the excitations. These 
functions can be constructed at both weak and strong coupling O IMl E31 EH ES] and given their 
knowledge in these regimes we will derive explicit expressions for the dispersion relations. Finally, 
note that the subleading, ©(log'' S), correction to the vacuum distribution density can also be casted 
into the form ()3.40p . as shown in [53]. The source terms in that case can be written as 

^ ^ 1^00 d,j^i2gt)Mm-gt5^, ^ ^ ^^^^^^ 

They differ slightly from the ones given in [53], due to a different way of introducing dj^{t). The 
relation between 'y±.{t),6^^{t), and the vacuum distribution density is given in Appendix [B] for 
completeness. 



Note that since the vacuum distribution density is even, there is no need for 7^(f), which are indeed zero due to 
=0. 
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To conclude this subsection, let us point out that the coefficients 72n-i l2n have a simple 
spin-chain interpretation. They are directly related to the spin-chain conserved charges. Indeed, it is 
not difficult to prove that 

(_l)n-l (-1)""^ 
<l2n = ^^-^2^l2n-l, ^L+l = 2^^^72. , (3.43) 

where the notation q* stands for the individual contribution of the excitation to the total charge 
(above the vacuum) qr, with qr = J2-kQr- The relations (j3.43|) are also valid for the vacuum solution. 
Namely, we have (g^+i = 0) 

and similarly for the vacuum subleading correction Sqf^. We recall that the spin-chain total charge 
Qr admits the decomposition 

Qr = qf (log S - + 5qf +Y.ql, (3.45) 

where the first two terms stand for the leading and subleading vacuum charge, previously denoted 
^vacuum -^^ (|3,l6p . Finally, generating functions for q* can be introduced as in (I3.30p . leading to the 
identities 

qt{2gt) = -2jl{2gt) , (t+C^gt) = -2^X{2gt) , (3.46) 

and similarly for qf and 6qf . The physical meaning, if any, of the functions ")X{t) and 7* (t) remains 
obscure, however. 



3.4 Fermionic subtleties 

The analysis for the fermions turns out to be subtle, as mentioned before, and requires additional 
comments. Fixing a fermion at a rapidity Uf, one could ask what are the coefficients k^, k^, that source 
the equations (|3.4U|) . Looking for them in Appendix ([B]), one would find two possible expressions, 
denoted kJ{, and k^^, R^^, for 'large fermion' and 'small fermion', respectively. The reason for this 
doubling is that the 'physical' rapidity for a fermion is the string spectral parameter Xf, not the 
spin-chain one Uf. Fixing Uf leaves two possible values for Xf, as it is well known. These are the two 
solutions to the quadratic equation 

Ui = Xi + g^ /xi. (3.47) 

The two kinematical domains to which we refer as large and small fermion in this paper correspond 
to Xf > g^ and < g^, respectively. In each of these domains, we can invert (|3.47p leading to 

Xi =\{uf + \Juf- (2^)2) , 

Xf=i(nf-y'nf-(25)2), (3.48) 

for large and small fermion, respectively. Of course, one would like to work directly in terms of Xf and 
forget about the separation between large and small rapidity domain. However, it turns out to be 
simpler to parameterize the dispersion relation in terms of the Bethe rapidity Uf; hence the necessity 
to deal with two set of coefficients, kJ{ and k^^, k^^. 

What happens at the level of the system of equations (j3.4Up when Xf crosses the boarder between 
the two domains is that the coefficients k'^, k'^ (with '-f = 'If or 'sf'), are no longer suppressed at large 
n. More precisely, they contain contributions of the type ((7/xf)", (xf/g)", for large, small, fermion, 
respectively. To get ride of these terms, one can redefine 7^(2(7t),7l':(2(7i), yielding a new system of 
equations, still of the type ()3.40p . but with a different set of coefficients k^, k^. This is the way one 
can switch from one domain to the other. In more sensible terms, the equations for large and small 
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fermion can be obtained from two different manners of looking at the scattering phase S4f. From the 
weak-coupling point of view, the domain > is the most natural one. In that case, one factors 
out the one-loop scattering phase and write 

= + f;-^:f 7' . (3.49) 

u — Uf — I 1 — g'^ /x^Xi 

Then proceeding along the lines presented before, one computes the inhomogeneous terms Xf(ii, Uf), 
Xf(n, Uf), by means of (j3.26p . and after a few of algebra get to the system of equations (|3.4U|) with 
source terms kJJ^, kJ^, as quoted in ([RM IR^ . But, obviously, the expansion around the one-loop 
scattering phase is not a good starting point for a fermion with x'j < g^, as this kinematical window 
shrinks to zero size at weak coupling. Instead of this, one would look for an expansion around Xf = 0, 
which means expanding around a fermion at rest. At exactly X{ = 0, the only effect of introducing a 
fermion is to reduce the spin-chain length by one unit P|, meaning that 

x+ 

54f(u,Xf = 0) = — . (3.50) 

X 

As it is the leading order contribution at small Xf, let us factorize it and rewrite 

54,f(n,Xf) = — fi^^") . (3.51) 

\ 1 — X{/X I 

Then, applying (|3.26|) , one finds that the inhomogeneous terms Xsf(n, Uf), Xsf(ti,tif), are simply given 
by 

f I — X{/x^\ 

Isf{u,Uf) +Is{{u,U{) = -i9«log — . (3.52) 

\l- Xi/X J 

In that form it is clear that at Xf = there is no deformation of the vacuum solution induced by the 
fermion, since 1,1 = 0, in agreement with the special role played by this point from the symmetry 
point of view [2]. Performing the steps explained before, starting with the expression (|3.52p . one 
would easily derive the set of equations (I3.40p with coefficients k^, ii^, as quoted in (|B.35|) . For 
consistency, one should observe that the two descriptions given above match when ~ ^f^. We shall 
see that this is indeed the case for the two expressions we will obtain in the large and small fermion 
domain for both the energy and momentum. 



4 Dispersion relations 

In the previous section, we derived a set of equations that control the contribution to the distribution 
density brought by the addition of an arbitrary excitation on top of the GKP string. In this section, 
we will explain how to construct a representation for the energy and momentum of the latter impurity. 

4.1 General expressions 

The GKP energy for an excitation with rapidity is given by the formula 

E^ = twisU + 2g'^q^, (4.1) 

where stands for the contribution of the excitation to the spin-chain energy, and twisty = 1 for 
•k = scalar, fermion, gauge field, while twisty = i for ★ = bound states of gauge fields {£ > 1). The 
identity (j4.ip can be converted into 

= twisty + 4(77* , (4.2) 

by using ()3.43p . So solving the system of equations ()3.40p for 7* would give us the energy, as expected. 
What about the GKP momentum? In principle, to obtain the momentum, we need to look at the 
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effective Bethe ansatz equations for the type of excitation we consider and extract the quantity that 
is conjugated to the length 21ogS. We did this analysis and derived the expressions given below for 
the momenta of the various excitations. Incidentally, we observed that the GKP momentum can also 
be obtained from the large u asymptotics of the odd part of the distribution density a{u). Namely, 
we found that 



at large u, with (5"j,(n) the correction to the total density sourced by the excitation. The identity (j4.3p 
can be rephrased in terms of the Fourier-like transform of o"j,(u), see Eq. (|3.32|) . as = 20(0). 
Now, using the relation between Cl^{t) and the auxiliary functions J±{t) (see Appendix [B|), one easily 
deduces that 

For any excitation, we have p°^'^~^°°P = 2u^,, except for a small fermion, p°'^'^~^°°^ (u^,) = 0. We conclude 
that to find the dispersion relation, we need to evaluate 7i,7i, and then apply (j4.2l 14. 4p . 

At this point we are still facing the difficulty to solve the system of equations (I3.40|) . Fortunately 
we are only interested in extracting the components 7i,7i, of the corresponding solution. For them 
it is possible to derive a representation in terms of 7±(t). The specificity of the latter functions 
is that they solve the system of equations (|3.40p with the inhomogeneous terms = 2g6n,i- As 
explained in [5B], this observation allows us to use the functions 7±(t) to project out the components 
7^,7j|', for functions 7^(t), 7j_(i), satisfying the same system of equations (|3.40p but with different 
inhomogeneous terms Going along these lines [3S], we find that 

457^ = E 2(2n - Ihl-A-i " E 2(2n)72^„4n , 

n>l n>l 

4573^ = E 2(2n - Ihfn-A-i + E 2(2n)72^n^2n ■ 

n>l n>l 

As we know explicitely the coefficients these identities provide a parameterization of the 

dispersion relations in terms of the solution to the BES equation for the large spin vacuum distribution 
of roots. 

For the cases we will consider, it is possible to find more tractable integral representations, starting 
from (j4.5p . Indeed, the coefficients are typically given by some integrals over the Bessel 

functions Jn(t). It suggests to commute summation and integration and then make use of (I3.39P to 
form the functions 7^(t)@ Simplifying a bit the obtained expressions, with help of identities for the 
Bessel's functions and using the equations for 7^ (t) when necessary [38] , one arrives at the dispersion 
relations given below. 

Scalar 

The dispersion relation for a hole carrying the rapidity u reads 

E{u) = l+ — -— (e / cos(nt) - 11 - / —7^(25^), 

Jo r e^-i V ^ Jo t (^ g) 



00 



dtj^{2gt)+j^{2gt) 



p(u) = 2u- , — ^ sin (ut) . 

^ ' Jo t e*-l ^ ^ 

Integrating the energy in ()4.6p with an even distribution density for the hole rapidities would lead to 
the representation of the generalized scaling function proposed in [38j. Taking the derivative of the 
momentum in ()4.6p with respect to the rapidity gives the leading-order contribution to the vacuum 
distribution density at large spin. 



-'vacuum 

vr 



(u) = ^-^p'iu) + 0(log° S) . (4.7) 



Commuting summation and integration is permitted by the fact that the Bessel's functions Neumann series H3.39P 
for "/^(t) are uniformely convergent over the real axis, because the coefficients 7^ are suppressed at large n. 
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This equality is due to the fact that the effective equations for the holes rapidities are obtained via 
the counting function Z{u) ^ p for the main roots, 

27rnh = Z{0) + 2tt dv p{v) = Z(0) + 2p{uh) log S + 0(log° S) , (4.8) 
Jo 

where we used that p{v) = Pvacuum(t') + 0{log^ S) and where is a hole mode number. The 
identity (j4.7p impHes that constructing the momentum of a scalar for arbitrary rapidity u is equivalent 
to solving the BES equation. 



Gauge field and bound states 

The dispersion relation for a gauge field (^ = 1), and its bound states {i > 1), carrying a rapidity u 
is given by 



+ 



E{u) 
p{u) = 2u 



dt 
T 
°° dt 
T 



1 -e-* 



e*-l 



(cos(nt)e-«/^-l 



7?(25t) ^ lf{2gt) 



1 -e- 



e* -1 



sin {ut) e-«/2 . 



(4.9) 



When u = 0, we have p{u = 0) = 0, and the expression for the mass E{u = 0) is in agreement with 
the findings of [20]@ 

Large fermion 

The dispersion relation for a large fermion with rapidity u is 

dtjf{2gt)-j^{2gt) 



E{u) = 1 + 



t e*-l 
+ cos {ut)^^{2gt) 



cos {ut) — 1 



t 



(4.10) 



/■°" dt 7g(2g^)+7g(2gt) 1 /--dt 



p{u) = 2u 

Note that this representation assumes that v? > (2g)^. At = {2g)'^ the two integrals 

/ cos(nt)7?(25t), / -sin(^xt)70(2gt), (4.11) 
Jo t Jo t 

have square-root branch points, pointing toward the fact that the physical rapidity is x = {u + 
\J'u^ — {2gY)/2, not u, as discussed already. The analytic continuation of (|4.1U|) in the complex 
It-plane will be explained later. 



Small fermion 

The dispersion relation for a small fermion with rapidity u is given by 



E{u) = l-- -^cos{ut)jf{2gt), 

^ -^0 * (4.12) 



1 dt 
PM = n/ -sin(nt)7?(2(?i), 

2 Jo t 



^®In the notations of [5D], the index £ is related to m, which counts the number of roots U3 in the stack £, by the 
equahty £ = m — 1, 
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assumme; that > {2gf. It can be written as a Taylor series around x = 0, which for x 
\{u — yju^ — (25)^) maps io u = OQ. It reads 



E{x) = l-Y.{-lT^l 



n>l 



p(^) = -E(-ir72^n.^i [^g 



n>l 



(4.13) 



with the coefficients 7^ defined by the expansion of 7^(t) over the Bessel's functions (|3.39p . In 
particular we have 

E{x) = I + 0{x^) , p(x) = + 0(x3) , (4.14) 



where we used that rcusp(5') = '^97 f ■ The asymptotics (|4.14p shows that the fermion has zero 
momentum at a; = and that it is mass 1 exactly to all loops, E(x = 0) = 1, in agreement with the 
discussion of [14J. The fact that the cusp anomalous dimension appears in the small x behavior of 
the fermion momentum (14.140 should not surprise us too much. Indeed, this is a direct consequence 
of the coupling between a fermion and the main roots, which comes through the spin-chain conserved 
charges. At small x, we have 

n ^ ~ = exp (-iP - iQ2X + 0(x^)) , (4.15) 

j=l ^ ~ ^4,j 

where P, Q2, are the total spin-chain momentum and energy. At large spin S ^ 1, the energy Q2 is 
dominated by the contribution of the vacuum distribution of main roots, 1g^Q2 = 2rcusp(5) log S" -|- 
©(log'^S'), such that the GKP momentum (which is conjugated to the length 21ogS') reads p(x) = 
^ cvlsx,{())x j (2g^) + 0(2;^). Expanding further in powers of x and using the relations (|3.44|) between 
vacuum conserved charges and coefficients 7^_i, one would easily reproduce the Taylor representa- 
tion (|4.13|) for the momentum. 

4.2 Weak coupling expansion 

At weak coupling, the dispersion relations (|4.6l 14. 9t I4.10[ I4.12|) can be evaluated by expanding the 
functions 7^(t) and 7?(t) over the Bessel's functions, see Eq. (I3.39p . We just need the expressions 
for the expansion coefficients 7^, see (|3.39p . at weak coupling. They are easily obtained in the form 
of an expansion in by solving iteratively the system of equations ()3.40p for 7^, i.e. with source 
terms (j3.4ip . Doing so, we find that the coefficients 7^ start as 



7f = 2g\l- + + 0(/) , 73^ = -^5' + OW) , 

Tf = - 4 ( y Cs + IOC5 ) 5' + 0(5') , 7f = + 0(5') , 



(4.16) 



and otherwise satisfy 7^ = 0{g"'^'^) for n > 1. Here Cn denotes the Riemann Zeta function ((z) 
evaluated at z = n. With the help of (I4.16p . one can compute the dispersion relations parametrically 
in terms of the rapidity u up to 0{g^) (i.e. three loops for energy and four loops for momentum) o 
The results are given below, where, in order to save space, we have introduced the notations 

V'i"^) {k, u) = l- {ipn {k + iu) +il)n{k- iu)) , 

• (4.17) 
V'i ^ (A;, «) = 2 ("^^ (^ + -i'nik- iu)) , 



^^Actually, the expression given in (|4.16p for the coefficient jf is not necessary for the evaluation of the dispersion 
relations at 0{g^). 
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with ^pn{z) = d"^^ logr(z). We also use 



(4.18) 



for the cusp anomalous dimension. 
Scalar 

At weak coupling, the dispersion relation for a scalar is given parametrically as 



E(u) 



p{u) 



1 + r 



2/ [4+) {lu)+6Cs 



+ 



a 



3 

2u + T. 



cusp 



{9) % 



(-) 



+ 



5 



^ k") (i, + 2vr24-^ (in) - 24C3Vi-^ (i ^) 1 + O {9') • 



(4.19) 



The expression quoted in ()4.19|) for the energy is in agreement with the three-loop computation of [36j 
using the long-range Baxter equation [37]. Moroever, the mass of a scalar E{u = 0) is directly related 
to the first contribution to the generalized scaling function fi{g) [171 I19j . Evaluating E{u = 0) 
with ()4.19p and simplifying the expression, one finds agreement with the weak coupling expansion of 
fi{g) [19J. This matching is actually guaranteed to all loops, by comparing the all-loop expression 
we found for E(u = 0) with the formula given in [38] for fi{g). Finally, note that it is possible to 
reexpress the momentum in terms of more familiar hyperbolic functions (at least for the first few 
terms in (I4.19P ). For instance, for the momentum of a hole at two- loop order, one has 



p{u) = 2u — 2TTg'^ tanh (vrn) + 0{g^ 



(4.20) 



Gauge field and bound states 

The dispersion relations for a gauge field {i = 1) and its bound states {i > 1) are given at weak 
coupling by 



Eiu) 



i + rcusp(5) [i'o'^ {s, u) - Vo(i)J - 2/ [vr' {s, y) + 6C3 

' ) {s, u) + 27r24+) {s, u) + 24C3^/'!+^ (s - 1, n) + 8 (^^^3 + SOCs) + O (/) , 



',(+) 



3 



p{u) = 2u + rcusp(5) 4 ^ (s - 1> - ' (s - 1, y) 



,4„/,{-) 



+ 



4"^ {s -l,u) + 27rV^^ {s - 1, u) - 24C3^i"^ (s, u) + O (g^) , 



(4.21) 



where s = 1 + i/2. For u = 0, we have p{u) = 0, and the masses E{u = 0) are in agreement with the 
results of [Snisij. 
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Large fermion 

The dispersion relation for a large fermion is 

E{u) = 1 + rcusp(5) Wo^^ (1, u) - M^)] - 2/ [4^^ (1, u) + 6C3' 



5' 
+ — 
3 



(1, n) + 27r''i;^^> (1, n) + 24(3^^' (1> ^) 



.2,/,(+) 



(+) 



= 2u + rcusp(5') V'o ^ (1. "w) 



cusp 



(5) 



2u 



12C3 



(^2^3 + 3OC5) 1 + O [g^) , 



+ 



4"^ (1, u)-^ + 27r^iP^-^ (1, n) + — - 24C3Vi"^ (1, u) 



+ 5' 



(4.22) 



The singularities at small u signal the presence of the square-root branch cut associated to the map 

^2 



hu + ^u^-{2gr) = u-^ + 0ig' 



(4.23) 



When expressed in terms of the rapidity x, the energy and momentum are regular for any real value 
of verified below. 



Small fermion 

The dispersion relation for a small fermion is 

E{u) = l + ^g^ + 0{g' 



p{u) 



2_l 2_l(^_f\ 



(4.24) 



Note that the weak coupling expansion is taken at fixed u. It means that the rapidity x scales as 

■,2 



x = i(n-v'n2-(25)2) = ^ + 0(/). 



(4.25) 



Here again the singular behavior at ti ~ indicates the emergence of the cut with branch points at 
u = i2g. But in distinction with a large fermion, we note that there seems to be no other singularities 
in the ti-plane of a small fermion. This is actually correct and both energy and momentum are 
holomorphic functions of x for |x| < g. 

Interpolating large/small fermion 

Let us verify here, at the level of the weak coupling expansion (j4.22lHT2^ . that large and small fermion 
are analytically related to one another. To this end, we start with the dispersion relation (j4.22p for 
a large fermion and reexpress it in terms of the physical rapidity x = u + g"^ /u. We find to two loops 
that 



E{x) = i+reu,p(5) 



4+)(l,x)-Vo(l) 



25^ 



4+)(l,x)+6C3--4~^(l,x) 



+ 0(/), 



p{x) = 2x + Tcuspig) V'o ^ (1) ^ 



(-)(l,x) + -f24+)(l,x)-^ 



+ 0(/). (4.26) 



In spite of their appearance, all the terms in square brackets are regular at small x, and keeping only 
the leading contributions, we get 



E{x) = 1 + (452^3 + 0(/))x2 + 0{x^) , p{x) 



rcusp(g) 

2g^ 



X + 0{x^ 



(4.27) 
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in agreement with ([12HS25I). We could have started from the small fermion domain as well. From 
the Taylor series for energy and momentum given in (|4.13|) and with the help of 




(4.28) 



one finds 



E{x) = l-4g^ 5](-l)\2n+ix2" + 0(/) 



l + 452(4+)(l,x)-V'o(l)) + 0(5') 



n>l 



p{x) = 2x + 4g^ ^(-l)"C2nx2"-i + 0(/) 



2x + Ag^4 \l,x)+0{g^) 



(4.29) 



n>l 



which is consistent with ()4.22U4.23[) . since rcusp(5') = 45^ + 0{g'^). 
4.3 Complex rapidity plane 

To go beyond the weak coupling expansion and improve our understanding of the complex rapidity 
plane, it is helpful to recall few facts about the solution to the BES equation and derive alternative 
representation for the dispersion relations. After these steps, we will show that fermions in the small 
and large rapidity domain are analytically related to one another. Then we will discuss what are the 
genuine singularities in the complex rapidity plane of the various dispersion relations and derive their 
large rapidity asymptotics. Finally, we explain why the higher-twist excitations can be thought as 
bound states of gauge fields by constructing a set of fusion identities for their dispersion relations. 

Alternative representation 

The prerequisite for the exploration of the complex rapidity plane is a brief reminder of the main 
properties of the solution to the BES equation \63\ [Ml [38l 165] . A remarkable account of the an- 
alytical properties of the functions 7±(t), which are directly relevant for our investigation, can be 
obtained from [M], where the solution to the BES equation was constructed in terms of resolvants. 
El Within our approach, most of these properties can be understood as being inherited from the 
Bessel's functions over which 7±(t) expand in a very controlled wayl^ The main point is that the 
Neumann series (|3.39p are absolutely convergent and converge uniformly toward the functions 7±(i) 
in the complex t-plane. It follows that the functions 7±(i) are entire functions over the complex 
t-plane, with an essential singularity at t = oo. They are also bounded over the real t-axis which 
allows us to perform the transformations 



converging for Im (n) ^ 0. Differentiating (j4.30p with respect to u and performing few rescalings give 
two of the resolvants of |64j . 

From the analysis of [54] or from the expansion over the Bessel's functions, we learn that A-i-(n) 
have square-root branch points at u = i2g, such that locally around these points 



with a±{u), f3±{u), being analytic at ti = i2g. The functions X±{u) can furthermore be continued 
to holomorphic functions of u in the complex plane C U {oo} outside the cut ^ (2<7)^- Explicit 

^^It is not difficult to establish the mapping between the resolvants of [B3] and the functions we manipulate in this 
paper. The dictionary is given in appendix of 64 and it can be applied here after performing an overall rescaling: 
7^(t) = 4(;7±(t) with 7±(t) the notations used in the dictionary of [64) . Similarly for the functions r^(t) (introduced 
below) and r±{t) in [Bl]. 

^®The expansion coefficients 7^, solving the BES equation in the form H3.40l[X^T1) . are well suppressed at large index 
n ^ 1 (taken at a fixed value of the coupling), which follows from the asymptotics of the kernel of the BES equation |61| . 




(4.30) 




(4.31) 
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expressions for the functions A±(n) in the lower-half plane then follow from (I4.30p and from the 
properties (A±(u))* = X±{u*), with * the complex conjugation!^ In particular, for > {2g)^, the 
functions X±{u) are real and can be written as 



\+{u)= -^%{2gt)cos{ut), \-{u)= -7^(25t)sin(ut), (4.32) 
Jo t Jo t 

as a consequence of (|4.30p and of the equalities (|A.5p for Bessel's functions. 

It is easily recognized that the functions A-i-(u) permit to analytically continued energy and mo- 
mentum of a small fermion. Indeed, comparing (I4.12|) and (I4.32|) . we conclude that 



1 



1 



(4.33) 



For later purpose, we also note that 



Esf{u) + Ps{{u) 



where 



X{u) = \+{u) - A_(u) 



jj^{2gt)e 



iut 



7^(t)^7?(t) + n^(i), 



(4.34) 



(4.35) 



where, in the first equality, the integral representation is for Im (n) ^ 0. One could as well consider 
the quantity Es{{u) — psi{u), which can be obtained from (j4.34p through the parity u — ^ — 
Now, applying what we have learnt before, we conclude that both energy and momentum of a small 
fermion are holomorphic functions of the rapidity u outside the cut ^ (^s)^- Equivalently, energy 
and momentum of fermion are holomorphic function of the Zhukowski rapidity x in the disk |x| ^ g. 

Around n = oo (i.e., x = 0) we get from (|4.3Up that X+{u) ~ and A_(u) ~ l/u, since 

7+(i) ~ t'^ and 7?(t) ~ t, implying that E^f = 1 -|- 0{p'^f) when ~ 0. When we go away from this 
point and get closer to the cut ^ (2^)^, the momentum typically increases in modulus, and, if done 
along the real «-line, the energy will also increase|fl Note finally that energy and momentum are 
generically complex around the cut, with values above and below the cut being complex conjugate of 
one another. So in order to increase energy and momentum further, one needs to enter the cut and 
change sheet. This is the connection with the large fermion domain. 

In order to verify that small and large fermion domain are correctly glued together, we need to 
analytically continue the functions X±{u) through the cut. This will be done later. For the time 
being, we simply report on the additional information we need to perform this step. This additional 
information comes from the fact that the functions 7±(t) solve the BES equation. A more useful 
statement is that they solve the equations 



°° dt 
T 

°° dt 
~t 



jfi2gt) j?{2gt) 



1 



e*-l 



7?(25t) , 7?(25t) 



1 



-1 



(cos {ut) - 1) = , 
sin {ut) = 2u , 



(4.36) 



that hold for < {2g)'^ . The system of equations (|4.36p is not equivalent to the original one (|3.4Up 
because it admits an infinite number of solutions. It is a sort of off-shell extension of the BES equation 
that turned out to be very useful in constructing its solution at strong coupling EH E5] • When 



^°We recall that the functions 7^(i), which are integrated in H4.30p . are real, if g is real. 
^^Note that \±{-u) = ±\±{u). 

^■^We do not know how to prove these naive expectations from general properties of the functions 7±(t), but we 
observe that they are correct at both weak and strong coupling. Note also that if we approach the cut starting form 
u = oo and running along the imaginary ti-line, then the momentum increases in modulus, but, since it is imaginary, 
the energy is decreased. 
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supplemented with appropriate analytical properties for 7± (t) the solution to (14.360 is of course unique 
and identical to the one obtained from (j3.40p . Here we just need to note that the equations (|4.38|) 
apply to our problem. 

If we aim at simplifying energy and momentum of gauge field and bound states, another function is 
worth considering. It is denoted T^{t) and is introduced by means of the transformation [63 | [Ml [65] 

r^{2gt) ^ (1 + icoth(l)) 7^(25*) , (4.37) 

with j^{t) defined in ()4.35p . In terms of this function, the equations ()4.36p can be combined into 

roo r "I 

/ r^(25t)e™*-*^/4-r^(-2(7t)e-™*+*'^/^ = iiV2 , (4.38) 

and it also assumes < (2^)^. The equation (|4.38p is far more simpler to analyze than the previous 
ones (|4.36p . especially at strong coupling. The price to pay is that the function T^{t) has more 
complicated analytical properties. It has for instance an infinite number of poles along the imaginary 
axis, due to the transformation (|4.37p . The quantity of interest here is not directly T^{t) but its 
Fourier transform, that we introduce as 

A{u) = I ^r^{2gt)e'^' . (4.39) 
J 27r 

It is defined for any real value fo u except at n = ±2g where it has (integrable) square-root sin- 
gularities, A(ti) ~ l/\/\u =F 2g\ . These two points split the analysis of A^(u) into the two domains 
^ (2(?)^. Here we will focus on the interval < (2g)'^. 

Starting from the segment < (2g)^, the function A(ii) can be analytically continued in both the 
upper- and lower-half u-plane. Indeed, after writing 

A{u)= r0(2<7i)e»*+ / r0(-25t)e-™*, (4.40) 

Jo 27r Jo 2tt 

we observe that for < {2g)'^ we can apply the equation (j4.38p and get both 



A(u) = -^e-/^ dtT^{2gt)e 



oo 



K[u) = - /"'dtr^(-25t)e-*"*+^'^/^ 
vr v2vr 



(4.41) 



From the two identities in (j4.4ip . we conclude that the function A(n) extends to a holomorphic 
function of u in both the upper- and lower-half u-plane. We recall however that A(m) has (integrable) 
square-root branch points at = (2g)^, such that the present analytical continuation is defined in 
the complex plane with a cut along > (2g)^. Notice finally that the function A(n) has essential 
singularities at both u = ibioo. Around these points, A{u) admits expansion in inverse powers of u 
starting as 



2\/2 _^j^/4 



A(n) = e- 

vr 



2u 



(4.42) 



for u in the upper- and lower-half plane, respectively. The large u expansion (|4.42p is merely asymp- 
totic, however, since T^(t) has poles in the complex t-plane. 

We can now derive an alternative representation for energy and momentum of gauge field and of its 
bound states. Namely, starting from the general expressions ()4.9p . we form the quantity Ei{u) +pi{u) 
reading 



Ei{u) +pe{u) = 2u + 2e 

dt 



+ J_ r -\e-'^/^T'^{2gt) (e*(-+-)*-l) + I ^ {-2gt) (e~^("- 
2\/2 Jo ^ L ^ ' ^ 



-it)t 



(4.43) 
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where e = ^/2. It then drastically simplifies when written in terms of A(u) that leads to the integral 
representation 

Ee{u) +pe{u) = ^ / ^ dv A{v) , (4.44) 

J u—ie 

with the contour of integration as depicted in Figure O A similar expression can be obtained for 
Ei{u) — p£{u), using the reflexion u — )■ —u, and reads 

ITT /■"+*'^ 

Ei{u) - pe{u) = — dv A{-v) . (4.45) 

J u—ie 

The two relations (|4.44l I4.45P conclude this discussion. We will now come back to fermion and prove 
the continuation between small and large rapidity domain. 




4g 

Figure 3: Two realizations for the contour of integration in the integral representation (j4.44p . associ- 
ated to a rapidity u with v?' ^ (25)^, respectively. Here k = i fixes the size of the bound state, i.e., 
the end points uzizik/2 = uzizii/2 of the contour of integration. The latter is chosen such as to avoid 
the cuts in the function A, running from ziz2g to ±00 and depicted by the two wavy lines. 



Interpolating small and large fermion 

To prove that small and large fermion interpolate each other, we will analytically continue through 
the cut the energy and momentum of a small fermion and show that they become identical to those of 
a large fermion ()4.10p . The path we will follow is depicted in Fig. [H We give the proof for the energy 
only, which can be easily adapted to the momentum. We only perform the continuation by going 
through the cut from above, but the same conclusion would be reached by analytically continuing 
from below. 

We start with the energy in the small rapidity domain, depicted by point A in Fig. [H Hence we 
can write 

E,i{u) = l-h+{u), (4.46) 

with A+(n) as given in (j4.30p . We can easily move in the upper half-plane and evaluate the energy 
Es{{u + iO~^) right above the cut when u'^ < {2g)'^ . We can go through the cut by writing the energy 
as 

1 dt 

Es{{u + iO+) = 1 + -\+{u-iO-) - / —-ff{2gt) cos (ut) = Eif{u - «0+) , (4.47) 

2 Jo t 

where we used the explicit form of X^{u). Observe that in writing (I4.47|) we have anticipated that the 
function we are getting to is the energy in the large rapidity domain, thus changing the subscript 'sf 
into 'If at the end of the equation. In particular, we stress that ()4.47p is different from Esi{u — iO"*"), 
which is the value of the energy evaluated below the cut in the small fermion domain. 
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Figure 4: Sketch of the complex w-planes for energy and momentum of a large (lower sheet) and small 
(upper sheet) fermion. Crosses indicate square-root branch points. Except for the central cut, the 
upper sheet is free from singularities. The lower sheet contains an infinite series of cuts which extends 
along the imaginary axis, with a spacing between two consecutive cuts equal to iti. The dashed line 
indicates the path followed during the analytic continuation from the small to large fermion domain. 



We know that the function X+{u) evaluated in (|4.47|) is analytic outside of the cut < {2g)'^, such 
that we need to consider the continuation of the integral in ()4.47p to get -E'if(u) away from the cut. 
To this end we note that the latter integral can be expressed as 



, —jXi^gt) COS {ut) = / — -t — (cos {ut) - 1) - / —jf{2gt) , (4.48) 

Jo t ^ Jo t e^— 1 Jo t ^ 

by using the first equation in (|4.36|) . which applies here since < {2g)'^. This identity is telling us 
that the left-hand side of (|4.48|) can be analytically continued into a holomorphic function of u in 
the strip (Im(n))^ < 1, since this is a property of the right-hand side of (|4.48|) . So we can easily 
move in this strip and get to the point B in Fig. [U which stands for a generic point in this strip. 
For real momentum, we need to choose a rapidity u at B that is real and away from the cut (i.e., 
> 2g). Then we can further simplify ()4.47p . since in this case A+(u) restricts to the real function 
given in (|4.32|) . Combining it with (|4.48|) . we get the expression for the energy of a large fermion 
quoted in M.lOp . 

We have thus shown the interpolation between small and large fermion. The two domains, and 
their respective sheets, are connected through their common cut along < {2g)'^ . This cut originates 
from A+(u) in ()4.47p . We learnt furthermore that this is the only cut in the strip (Im (u))^ < 1 on the 
large fermion sheet. As we shall see, outside this strip, the energy of a large fermion has an infinite 
number of cuts, in distinction to what happens in the small fermion domain. 



Singularities in the complex rapidity plane 

We now discuss the singularities of energy and momentum in the complex u-plane for all the exci- 
tations. A natural starting point is to inspect the weak coupling expressions quoted in (j4.191 14. 2H 
14.221 I4.24p . We observe that both energy and momentum have poles along the imaginary u-axis, 
whose precise locations depend on the flavor of the excitation. For scalar, for instance, we find poles 
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at Up = i{2Z + l)/2, for both energy and momentum|fl We notice then that the orders of the poles 
increase with the loop order. This is suggesting that the weak coupling expansion breaks down around 
these points and thus does not give a proper account of the singularities in the u-plane (at any given 
order in g'^). 

It is then not difficult to figure out what the genuine singularities are: at any finite value of the 
coupling constant we only have cuts of size connecting square-root branch points. The reason 
being that we can decompose energy and momentum into sums of terms of the type 

X±{u + im), (4.49) 

with m belonging to some subset of Z or (2Z + l)/2, depending on the flavor of the excitation 1^ 
The weak coupling expansion of (14.490 produces poles of various orders at Up = —im, but, at finite 
coupling, the actual singularities of (j4.49p are two branch points at Ub-p = Up i: 2g, respectively. 
Summing over m leads to series of cuts in the complex -u-plane. From this perspective, the weak 
coupling expansion, taken at a rapidity u Up, can be thought as Laurent expanding around these 
cuts, thus generating poles sitting at their individual centers of mass. 




Figure 5: Complex li-plane for energy and momentum of scalar, gauge field, and bound states. Crosses 
indicate square-root branch points. Wavy lines indicate a choice of cuts that is natural from the weak 
coupling perspective. The dashed lines represent the real and imaginary u-axis, respectively. The 
integer k depends on the excitation. It is equal to 1 for a scalar and to £ for a bound state of i gauge 
fields. For a large fermion, the picture with k = 2 can be applied if completed with an extra cut 
along the real u-axis. For a small fermion, we would only keep the latter cut and remove the two 
semi-infinite towers in the upper- and lower-half plane. For all excitations, within the weak coupling 
expansion, the sequence of cuts degenerates into a sequence of poles along the imaginaris it-axis. 

Now, to complete our discussion of the singularities in the complex rapidity plane, it remains to 
list the center-of-mass positions Up of the cuts for all the excitations. They are the same for both 
energy and momentum of a given excitation. We find that Up = i(2Z + l)/2 for a scalar, Up = iL for 

^^For other excitations, it is interesting to note that a few new poles shows up as we increase the loop order. At 
0(cp^ however the pattern of poles is complete and found to be the same for both energy and momentum of a given 
excitation. 

^*In practice the sums are taken over an infinite subset of 1 or (2Z + l)/2, except for a small fermion where only m = 
appears. Then, in order to make these sums convergent, one should consider terms of the type \±(u + im) — A±(im) 
instead of H4.49[) . This does not affect the argument we are giving here. 
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a large fermion, Up = for a small fermion, and finally Up = zizi{i/2 + N) for gauge field {i = 1) and 
bound states {£ > 1). They are depicted in Fig. IHandEl As we have already mentioned, some of the 
expressions we are manipulating are directly related to the resolvants of [6l]- For them the pattern 
of cuts is in agreement with the findings of [SJj and otherwise stands for a quite straighforward 
extension. Finally, we recall that the central cut (at Up = 0) of large and small fermion are glued 
together and that it disappears when energy and momentum are expressed in terms of x. It means 
that all the singularities in the x-plane lie in the large fermion domain, |x| > g, and more precisely 
at Xb-p = x(ub-p) with Ub-p = ±2g + i{Z — {0}). 

Quasiclassical limit 

There is another singular point for energy and momentum in the complex rapidity plane, namely the 
point at infinity. There energy and momentum admit asymptotic expansion, whose leading terms 
turn out to be controlled by the few observables that parameterize the vacuum energy, i.e. the cusp 
anomalous dimension and virtual scaling function. 

At large rapidity u, with u going to infinity without crossing the cuts, energy and momentum are 
found to scale as 

E,{u) = T,^,p{g) (logn - ^(1)) + twist, + + 0{l/u) , 

2 (4.50) 

p,(n) = 2u- - rcusp(9) + 0{l/u) , 

where twist* = 1 for * = hole, fermion, gauge field, and twist, = i for bound states {i > 1). The 
results (|4.5U|) are given here for n ~ oo with Re (u) > and are easily extended to the case Re (u) < 
by using the parity of energy an momentum under u — >■ —u. Note also that in the case of fermion, the 
asymptotics (|4.50p holds for a large fermion, since for a small fermion the point u = cxd is associated 
to X = 0, that is a fermion at rest. Say differently, the asymptotics (I4.50p is a large momentum 
expansion, 

E^{p) = reusp(5) (log (p/2) - V(l)) + twist, + + 0(l/p) . (4.51) 

The logarithmic scaling of the energy (j4.50|) is a longstanding result [TSl Ej • The fact that 
it translates into the scaling ()4.5ip is a more recent observation derived at strong coupling and 
conjectured to hold exactly in [15j. It is directly related to the change of asymptotics [35 t [T6 l [T7] 

A-Sr. A2{g) logs ^A3{g) logs, A^ig) = SF^^M , (4.52) 

accompanying the entering of a hole rapidity in the quasiclassical domain u ~ S" ^ ijffi This 
large hole brings us far away from the vacuum and closer to some highly-excited trajectories in the 
spectrum of large spin anomalous dimensions. The string theory dual of these trajectories are the 
spiky strings [58], that generalize the GKP string to configurations with more than two fold-points. 
Accordingly, at strong coupling, the asymptotics (j4.5Up is observed in the giant hole regime jl5j . 
where the excitation is semiclassical and connected at large momentum to the spiky-string solutions. 

One can easily verify the asymptotics (|4.5Up at weak coupling by using the formulae (|4.19l 14. 2H 
I4.22|) , and given the weak coupling expansion for the virtual scaling function [19^ [53] 

B2{g) = -24C(3)/ + y (7r2c(3) + 30C(5)) / + 0(/) . (4.53) 

It is not difficult then, following the hint of the weak coupling expansion, to rewrite appropriately the 
integrands in the integral representations (|4.6l 14. 9t l4.1Up such that the result (|4.5Up becomes apparent 
at any value of the coupling. 

^^This regime is quasiclassical from the point of view of the Baxter equation [351 117| . 
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4.4 Prom gauge field to bound states 



We are now in position to prove the fusion identities that Unk energy and momentum of bound states 
to those of the gauge field. We recall that these relations read 



Egi{u + ji 

j=0 



i(£-l)/2) 



(4.54) 



with -E'gf(tt) = Ei=i{u), and similarly for the momenta. They are saying that the twist-£ excitation 
with index £ is a ^-string of gauge field. These identities are actually quite manifest when looking at 
the contour integral representations ([Hll [135]) , such that the main point is to explain why they are 
obscured at weak coupling. 
Let us focus on the simplest case of a bound state of two gauge fields. We have both £ = 2 and 



£^first-bs(^i) = Egf{u + I) + Egf{u 



(4.55) 



where £'first-bs(^i) = Ei=2iu)- The equality (14.550 immediately follows from the concatenation of the 
contours of integration when using (|4.44l I4.45|) to represent the right-hand side of (j4.55p . The same 
holds true for the momenta. On the other hand, one can easily verify that the identity ()4.55p is not 
satisfied by the weak coupling expressions (j4.9|) . The reason for this apparent disagreement is that 
the fusion relation ()4.55p assumes that the rapidity u lies in the strip —2g < Re (u) < 2g, while the 
weak coupling expansion works outside of it. Performing the shifts in (|4.55p inside and outside this 
strip are two non-equivalent operations. This is due to the presence of cuts connecting the branch 
points at M = ±2(7 it i(2Z + l)/2, in the complex rapidity plane for energy (and momentum) of a 
gauge field. 

To verify the previous assertion we shall perform the shift in (I4.55|) for arbitrary u, assumed real 
for simplicity. We get that 



Egf{u + I) + Egf{u - I) = -Efirst-bs(n) + -Rdisc(li) , 



(4.56) 



where the remainder function i?disc(^) is related to the discontinuity of A{v) through its cut, see 
FiglH An explicit expression for the latter function is 



-Rdisc(^i) 



°° dt 

T 



jf{2gt) 7?(25t) 



1 



-1 



(cos {ut) — 1) 



(4.57) 



By means of the representation (I4.57p . we verify that iidisc(^^) = when u"^ < {2g)'^, in virtue of 
the first equation in (j4.36p and in agreement with our previous discussion. However, the remainder 
function i?disc(^) is non-zero for u'^ > {2g)'^, and this is why the identity (j4.55p does not hold at weak 
coupling. A similar analysis carries over for the momentum, for which one can find 



Pgf(n + |) +Pgf(u - §) = PRmt-hsiu) + R, 



''disc(''^) ' 



where 



R'disciu) =2u- — 
Jo t 



7^(2gt) 2f{2gt) 



-1 



sin (ut) 



(4.58) 



(4.59) 



To some extent, the fusion identities ()4.54p still make sense in the domain (Re(n))^ > {2g)'^, but to 
interpret them correctly we need to specify that the gauge-field rapidities in (j4.54p are not all lying 
on the same sheet. Only one of them can actually be on the physical sheet, that is on the sheet that 
contains the real-momentum line. 
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Figure 6: Fusion of two gauge- field excitations into a bound state. For a rapidity u'^ < {2g)'^, depicted 
on tlie left, the two shifted contours of integration add up to form the contour of integration for the 
first bound state. For a rapidity > {2g)'^, depicted on the right, the two shifted contours of 
integration hit the cut and extract an extra contribution fom the discontinuity. 

5 Strong coupling regimes 

In this section, we analyze the dispersion relations at strong coupling. We will show that the latters 
admit three qualitatively distinct regimes, that parallel the ones found for a magnon over the BMN 
vacuum \27\ Wl\ I46j . There are the perturbative regime \13\ [H] with E,p ~ g^, the near- flat space 
regime with E,p ^ 9^^^-, and the giant hole regime [15] with E,p g. The near-flat space regime 
interpolates between the perturbative and semiclassical domain. The expansion parameter is 1/g in 
the latter two cases and in the former. 

In the perturbative regime and to leading order at strong coupling, the dispersion relations are all 
relativistic and the pattern of masses agrees with the one obtained directly in string theory [131 jffl 
In the near-flat space regime, all excitations are found to be massless with dispersion relation E = p, 
to leading order at strong coupling and at fixed E = E/ g^/^ ,p = p/ g^/'^. The leading-order dispersion 
relation in the giant hole regime is also independent on the flavor of the excitation. Its expression, 
derived below, agrees with the one obtained in [TS] using the finite gap technology. 

For scalar only, there is an extra regime of exponentially small energy and momentum, ~ e"'^^. 
The presence of this regime was predicted and elucidated in [13] and it originates from the non-trivial 
infrared dynamics of the transverse fluctuations in S^, which turns out to be controlled by the non- 
linear 0(6) sigma model. At the perturbative level, one has 5 massless Goldstone bosons associated 
to the spontaneous symmetry breaking of 0(6) down to 0(5). But the symmetry is restored non- 
perturbatively and one ends up with a 0(6) vector multiplet of massive particles with mass m ~ e~'^^. 
This non-perturbative regime was previously observed to emerge from the Bethe ansatz equations 
in [39t IBS} BO] in the study of the generalized scaling function [TO] . 

5.1 Scalar 

We start our analysis with the dispersion relation for a scalar in the non-perturbative regime, E,p ^ 
e~^^ . It appears when the Bethe rapidity u is kept of order 0(1) while taking g ^ 1. When 
the rapidity gets larger, but remains in the domain 1 <^ u <^ 2g (for right mover), we enter the 
perturbative regime with E,p 1. We will analyze it up to 0{\/g) correction. For n — 2/7 ~ 1, we 
fall in the near-flat space regime with E,p g^^'^- The perturbative expansion is reorganized in this 
regime and the expansion parameter becomes 1 / y^. We will construct the dispersion relation in the 
near-flat space regime up to the next-to-leading order. Eventually, for a rapidity u » 2g, we find a 

^®Up to the missing mass 2 boson. 
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semiclassical regime where E,p ^ g. This is the giant hole regime [15], that we wih study to leading 
order at strong coupling. 



Non-perturbative regime 

To deal conveniently with the non-perturbative regime, we extend and apply the representation 
of \3S\ providing energy and momentum in the form of hyperbolic Fourier seriesQ It reads 



E 



P 



E 

n>0 

E 

n>0 



(-l)"m4n+i cosh ((4n + 1)6) - (-l)"m4n+3 cosh ((4n + 3)6) 
(-l)'"m4n+isinh((4n + l)e) + (-l)"m4n+3 sinh ((4n + 3)0) 



(5.1) 



with 6 given in terms of the Bethe rapidity as 6 = ttu/2. The representation (j5.ip is exact but it 
assumes 6^ < (irg)'^, i.e. < (2^)^, for the sums to have a chance to converge. The coefficients 
m4n+i,m4n+3, are functions of the coupling constant only and they can be expressed in terms of the 
solution to the BES equation, see Eqs. ()D.19l [D.2ip . For the time being, we only need to note that 
these coefficients are exponentially suppressed with the coupling as rUn ~ exp {—mrg). 

It immediately follows that, at strong coupling and for a fixed rapidity 6, the dispersion relation 
becomes relativistic |14^ [38] 



E' 



p' = m\ , (5-2) 

up to subleading ~ exp (— 47r5) ~ m| corrections. The coefficient mi, entering ()5.2p . is known to 
admit the strong coupling expansion 



mi 



-■Kg 



3-61og2 

327r5 ^ ' ' 



k = 23/4vrV4/r(5/4) , 



(5.3) 



where the (embedding) parameter k was first obtained in [14j from direct string theory computa- 
tion [66] . The relation ()5.2p allows us to identify the six scalars of the gauge theory with the vector 
multiplet of excitations of the non-linear 0(6) sigma model, which stands for the low-energy effective 
theory of the GKP string [14j. 

To go beyond the leading order expression (15. 2p . it is convenient to first introduce the rescaled 
variables e = E/mi, p = p/nii. Then the dispersion relation can be found as an expansion in powers 
of mf at fixed p. For the first correction to ()5.2|) . we need to consider 



e = cosh 6 — cmf cosh (36) + 0{m'l) , 
p = sinh0 -|- cm^ sinh (30) + 0{mf) , 



where c = m^/mi. The relation between p and 6 is easily inverted as 



6 = arcsinhp — c 



3p + 4p3 

7i+7 



leading to the dispersion relation 



l + p2 



1+- 



+ 



ml + 0{m\] 



-m\ + 0{mi 



The coefficient c is evaluated using (102111 and reads 



c = 



ma 



m 



64(12)1/4^2^ 



+ 0(1/5' 



■^^This representation will not be deduced here from the general formula 
for the methodology. 



(5.4) 

(5.5) 
(5.6) 
(5.7) 

and the reader is referred to [381 165| 
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Looking at the equality (15. 6p we note the emergence of singularities at p = ±i. It is a fake actually, 
and it stands for a mass renormalization. In our units, the mass is equal to 1 at leading order, but it 
receives a correction ~ cmf . Taking this effect into account leads to 



icp^mf + 0{mf] 



(5.8) 



with the (rescaled) mass m = 1 — cmf + 0(mf). Now the factor in square bracket includes the genuine 
corrections only. The dependence on the momentum p of the first correction in (|5.8|) suggests that 
it originates from higher-derivative operators in the effective (low-energy) Lagrangian density for the 
0(6) excitations. These irrelevant interactions break the boost invariance in agreement with the fact 
that it is a low-energy symmetry only [14]. 

It is not complicated, though tedious, to include subleading corrections to (|5.8|) of order 0(mf") 
with n > 2. Doing so, one would observe that higher powers of p are generated at each steps. There 
is a bound, however, in the sense that the large p asymptotics of the 0(rn-^") contribution cannot 
exceed p^". It points toward the existence of a non-trivial scaling limit when p ^ 1 with prrii kept 
fixed. This is, of course, the string perturbative regime, with p = pnii = 0(1) and mi — )• 0, that we 
analyze below. As a final remark, we notice that in terms of the rapidity 9 the perturbative regime 
requires the scaling 

9 — ng + jlogg = fLnite , (5.9) 

as follows from (|5.5p and (|5.3p . for right mover (^ > 0). It means that 1 ^ <C irg, or equivalently 
1 <C n <C 25. 



Perturbative regime 



To deal with the string perturbative regime, it is convenient to first reparameterize the dispersion 
relation (|5.ip . keeping only the terms that are relevant when mi ^ p ~ 1. For right mover, it amounts 
to neglecting all contributions to energy and momentum which are powers of e~^ in ()5.ip . Doing so 
and introducing a new variable, z, as 



log mi -|- log z , 



(5.10) 



which is kept finite when g ^ 1, we get from 



P 



n>0 
n>0 



4n+l 



+ i-irdnz 



4n+3 



(5.11) 



with the coefficients, see Eq. ()D.2ip . 

4n+l 
1 



dn 



m 



4n+3 



-n-1 



m 



(5.12) 



Note that, by definition, cq = 1. The relations (|5.1ip are valid up to non-perturbative contributions 
~ m\/p, but they include all perturbative corrections in \/g. Strong coupling expansion for the 
coefficients Cn,dn, can be computed order by order in 1/g from the solution to the BES equation, 
and their leading expressions are given in ()D.2ip . Finally, we note that, at a given order in 1/g, 
the infinite sums in (jS.lip truncate, see Eq. ()5.12p . and the dispersion relation E = E{p) is easily 
obtained. 
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Applying (jS.lip up to 0{\/g) contribution, which requires expanding (jS.lip up to 0{z^), we find 
the dispersion relatiorl^ 

r -r^,^ \\ ^ 

p 



E 



P 



1 



8(12)5/4vr2g 



+ 0(1/5^) 



(5.13) 



The leading term, E = p, is in agreement with the string theory result |13j . Moreover, if interpreted 
as a low-momentum expansion, the relation E = p + 0{p'^) holds true to all loops, that is up to 
non-pertubative corrections ~ m\/p. This has to be attributed to the identification of a scalar with 
a Goldstone fluctuations in [Hj. 

For illustrative purposes, let us include the 0{l/g'^) correction to (j5.13p . They are obtained by 
working at 0{z^) in (jS.lip . The dispersion relation then reads 



E 



p 



9 



9 



(5.14) 



where the explicit values of a, ... ,6, are not relevant here. We see that a new pattern appears. 
Namely, increasing the order in 1/g comes with an increment of the maximal power of p by 4 units. 
This is not apparent if only the 0{l/g) correction is considered. It suggests a non-trivial scaling when 
p gets as large as g^^^, with the energy given by 



E -p = g' 



-1/4 



a 



P 

.1/4 



+ 5 



P 

jl/4 



+ ... 



+ Oil/g 



-3/4N 



(5.15) 



As we shall see, the regime p ~ g^^^ interpolates between the perturbative and semiclassical domain, 
and it is therefore the analogue of the near- flat space regime for magnon above the BMN vacuum [U] . 
The difference with the BMN case is that energy and momentum scale like g^^^ and that the expansion 
runs in inverse powers of 1/y/g. 

Notice that assuming that p ~ g^^'^, or equivalently z ~ g^^'^, means that 6 — irg is kept finite for 
g ^ 1. In terms of the Bethe rapidity it means that u — 2g = 0(1). It is important to observe that 
u gets not rescaled when working in the near-flat space regime. Indeed, we recall that energy and 
momentum for a scalar have branch points at u = 2g ziz i{2n + l)/2, with n E N (and similarly for 
u < 0). Rescaling u by a power of g would cause the cuts to collide, leading to singular behavior at 
u = 2g |64j . This is not what we want as we would like to transit between u < 2g and u > 2g. Later 
we will perform the rescaling u ^ u = u/{2g) in the giant hole regime u > 1. 



Near-flat space regime 

In the near flat space limit, we define ^ = 6 — irg that is kept fixed for ^ ^ 1. In terms of this new 
variable, the overlap with the large momentum limit in the perturbative regime is at ^ = — oo. It is 
also convenient to introduce rescaled energy and momentum as £' = E/g^^^ and p = p/g^^^. The 
dispersion relation can then be written parametrically as 

E = Xo{0-Xi{C)/V9, P = Xo{0 + Xi{0/V9, (5.16) 
where the functions Xo,iiO given by 

-. oo 1/4 oo 

Xo(0 = E(-l)""^4n+ie(^"+i)(-^+«) , MO = V E(-l)""^4n+3e(«)('^^-*-«) , (5.17) 

Incidentally, the expression H5.13P could have been obtained immediately by taking the large p asymptotics of (|5.8p . 
It does not mean however that there is a one-to-one correspondence between 0{l/g") correction in (|5.13p and 0{mi") 
contribution in (|5.8p . For instance, to get the full answer at 0{l/g^), from the asymptotics of (|5.8p . one would have to 
include corrections up to 0(mf). 
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with the coefficients m4„_(_i, m4„+3, the same as before. Note that the dispersion relation (15.160 is just 
a rewriting of (|5.11|) and thus is also exact up to non-perturbative corrections. The main difference 
with the perturbative regime is that we cannot truncate the sums in ()5.17p . since all contributions 
are of the same order in g at given ^, see ()D.2ip . The functions Xo,i(0 defined in (j5.17|) are both of 
order 0(1) at strong coupling fixed), and admit an expansion in 1/g. To leading order, we find, 

see dEMHEHEil]), 

^»'«' = * (' + "-'''-'V + 0(1/9) . 



These functions are positive and increase with They are exponentially small at large negative ^ 
and scale like some fractional powers of ^ at large positive ^. 

Since the functions Xo,i(0 ^-re of the same order in the coupling at given ^, it immediately follows 
from (|5.16|) that the dispersion relation in the near-fiat space regime can be written as 

E = p + 0{1/^). (5.19) 

This result is actually universal, in the sense that its holds true for any excitation (scalar, fermion, 
etc) in the regime ~ p ~ g^^^- The subleading 0(1/ ^/g) corrections are however different in each 
case. For the scalar, this correction is controlled by the functions (j5.18p . 

E-p = -2xi{0/^ + 0{l/g), P = Xo{0 + 0{l/V9)- (5.20) 

It does not seem possible to solve this parametric representation for generic ^. But it can be done 
easily for the two extremal values ^ = ±00. At large negative ^, see Eq. dEM]), the near-flat-space 
momentum p is small, p <C 1, and the 0(1/ y/g) correction to the near-flat-space energy scales like 
~ p^. This is just saying that the dispersion relation takes the form (|5.15|) . which is expected as 
^ ~ —00 overlaps with the large momentum limit 1 <C p (<C g^^^) taken in the perturbative regime. 
At large positive ^, using asymptotics quoted in ()D.27p , we flnd that the near- flat-space momentum 
is large and that the 0{l/y/g) correction to the energy E reads 



E-p 



1 /3p\^/3 

■5 T + 



/V9 + 0(l/s), p»l. (5.21) 



We note that for = — irg ^ 1, the Bethe rapidity u = 29 /tt has become larger than 2g. It 
means that we are entering the giant hole regime, which is the last destination of our trip. If the 
near-flat space regime does indeed interpolate between perturbative and semiclassical regime, the 
large momentum asymptotics (|5.2ip should match the small momentum expansion in the giant hole 
regime. In the latter domain, energy and momentum are of order 0{g) and it is better to work with 
£ = E/{2g) = E/{2g^/^) and V = p/{2g) = p/{2g^/^). In terms of these variables the identity (fOTI) 
can be rephrased as 

£ = V-^ {3Vf^^ + ... , (5.22) 

and it should have the meaning of a low momentum expansion, V <^ 1. This is effectively the correct 
expression [15], as proved below. 



Giant hole regime 

To obtain the dispersion relation for a giant hole, it is actually simpler to start with the formula ()4.6p 
and immediately take the limit g ^ 00 with u = u/{2g) kept fixed and > 1. Then, to leading 



41 



order, we need to evaluate 

dt 



£{u) ^ E{u)/{2g) = / (7?(t) - 7?(i) (cos {ut) - 1) + 0(1/5) , 
V{u) = p{u)/{2g) =2u-J^ ^ (7?(t) + l^{t)) sin {ut) + 0(1/5) , 

with V? > 1. These expressions are actually common to any excitation carrying rapidity > 1, 
meaning that the dispersion relation in the giant hole regime is flavor independent. Note that in the 
case of fermion, the formulae (|5.23|) are obtained in the large rapidity domain. 

The integrals (j5.23|) are most easily taken after differentiating with respect to u. The two constants 
of integration can be fixed afterwards with the help of the large u asymptotics (I4.50p and of the 
known strong coupling expression for both the cusp anomalous dimension [12j and the virtual scaling 
function [53j . Using the leading order solution to the BES equation [62], 

lf{t)+^^^{t) = f^dv (^l±^y^'e-»* +0(1/5) , (5.24) 
one easily obtains that 

where we have dropped the 0(1/^) loop corrections. Now to integrate these expressions, it is helpful 
to switch to the Zhukowsky rapidity x = u + ^/v? — 1. One gets 

8 = j^^^l-l/z^ + Ci, V = 1^ dz^l-l/z^ + C2. (5.26) 

The correct choice for the constants of integration will turn out to be Oi = 02 = 0. Setting them to 
zero and performing the integrals, we obtain the semiclassical dispersion relation as 




1 - l/x2, 

(5.27) 



- arctan^xy^l — 1/x 

A simple reparameterization leads to the dispersion relation for a giant hole found in |15| . 

The dispersion relation (j5.27|) implies that both energy and momentum are small for x ~ 1"*". 
Expanding around this point, we verify that at low-momentum the energy scales as |15j 

£ = r - ^{srf/^ + o{r'^/^) , (5.28) 

confirming the overlap with the large-momentum limit in the near-flat space regime, see Eqs. (|5.2H 
I5.22|) . In the opposite limit, x — )• oo, energy and momentum are large. They read 



£ = log{Au)-l + 0{l/u^), P = 2u-| + 0(l/u), (5.29) 

where we have restored the (rescaled) Bethe rapidity u = x/2 + 0(l/x), with u ^ 1. Rescaling 
energy, momentum and rapidity, the identities (I5.29P become 

E = 2g{log{2u/g)-l) + 0{l/u'^), p = 2u - irg + 0{l/u) . (5.30) 

This is in agreement with the general expression ()4.50p for the large u asymptotics of the dispersion 
relation, given that the cusp anomalous dimension and virtual scaling function, entering (j4.50p , are 
given by [13 [13 

rcusp(5) =25 + 0(l), 52(5) =45 (log (2/5)+ V'(l))- 45 + 0(1). (5.31) 

We verify, by this way, that the choice Oi = O2 = for the constants of integration in ()5.26p was the 
correct one. 
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5.2 Fermion 



Our investigation of the dispersion relation for a fermion will be mostly restricted to the perturbative 
regime. This is the regime that is connected to the point x = (fermion at rest) and it is found 
in the small rapidity domain. It means that we shall start with the formulae ()4.12p . The strong 
coupling expansion is taken at a fixed value of the rescaled rapidity u = u/{2g), with v? > 1 
for a real momentum fermion. For a small fermion it means that the physical (rescaled) rapidity 
X = x/g = u — ^/v? — 1 fulfills < 1. Under the same assumption on a large fermion would have 
rapidity x = u + ^v? — 1 satisfying xP' > \. This is the giant hole regime, with the same dispersion 
relation as found for a scalar, to leading order at strong couplinglfj The near-fiat space regime 
controls the transition between xP' < 1 and xP > 1. Its analysis will not be done thoroughly, but we 
will verify that the leading-order dispersion relation E = p + 0{1/^) holds true with E = E/g^^^ 
and p = p/g^^^ kept fixed at strong coupling. 



Perturbative regime 



For evaluating energy and momentum of a fermion in the perturbative regime, we start with the 
representation (j4.34p and look for a strong coupling expansion, assuming u = u/{2g) fixed {v? > 1). 
It is obtained by expanding at large g and fixed t the function 7^(t). Using the known strong 



coupling expansion for the latter function 
the dispersion relation in the form 



and the formula (|4.34|) . one easily obtains 



E- 



■p 



1 



fix) 



with the function f{x) given up to 0{l/g'^) in Appendix [Dl Since {E,p) 
transformation x — )• —x, energy and momentum read 



E 



p 



1 + x 

1 - x 

1 + x 

1 - X 



1/2 



fix) + 



1 



1/2 



fix) 



1 + X 
1 - X 

1 + x 



1/2 



/(- 



1/2 



/(- 



-X] 



-X] 



(5.32) 

{E, —p) under the 
(5.33) 



while the dispersion relation is most easily obtained by means of 



E' 



= fix)fi-x) 



(5.34) 



Notice that the representation ()5.33p should not be interpreted as saying that energy and momentum 
have branch points at a;^ = 1. This is actually the case to any order in the strong coupling expansion, 
but, as we already said, energy and momentum are regular for x real. What happens is that the 
strong coupling expansion is singular at a;^ = 1. All the branch points, lying in the large fermion 
domain, collide at strong coupling in the rescaled rapidity x = x/g. They accumulate at = 1 
generating the illusion that energy and momentum are singular at = 1. The resummation of the 
1/g expansion would reveal that this is not the case [64]. 

After this brief remark, let us come back to the dispersion relation (I5.34p . To leading order at 
strong coupling, one easily find that f{x) = 1 + 0{l/g), giving the relativistic dispersion law of a 
mass 1 particle, 

E^-p^ = l, (5.35) 

in agreement with string theory [13]. The relation between momentum p and rapidity x takes the 
form 

p 



^^l+P^ 



1 + 0(1/5) 



(5.36) 



^® Indeed, starting from the representation (I4.10[) for the dispersion relation of a large fermion, and taking the strong 
coupling limit at u fixed (-u^ > 1), one would get the same expression as the one derived in (|5.23|) for a scalar with 
rapidity -it^ > 1. 
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meaning that x is the relativistic velocity of the fermion. It is not difficult to include higher-loop 
corrections. To 0{l/g'^) accuracy, one finds the dispersion relation (see Appendix [D]) 



with 



(5.37) 



K 31og2 1 4K 1 4K 



where K is the Catalan's constant. The first correction to the {x,p) relation reads 

P 



with 



1 + (a + bp')/g + 0{l/g' 



(5.39) 



3 log 2 1 , 31og2 

Looking at (|5.37|) . we verify the absence of correction to the mass, E{p = 0) = 1. We also observe 
that the identity E{p = ±1) = is satisfied [14] . It is so to all orders at strong coupling. We note 
that we can reach the points p = iti by sending x to infinity along the imaginary axis, x = itioo, 
see Eq. (|5.36|) . To get to these points we can think starting at x = 0, at which p = and E = 1, 
and then move along the imaginary axis. At intermediate values x = itz, we are at p = ibi/-v/2 and 
the energy has decreased up to £^ = l/\/2- At these points, we are about to leave the small fermion 
domain, since in terms of the Bethe rapidity we are right above or below the cut at u = ±0. If we go 
further we enter the large fermion domain and thus taking x — )• itioo requires analytically continuing 
through the cuts in the upper- and lower-half plane of this domain. Here the cuts are not apparent 
since, as we have said before, they have collided in the current scaling. It should be clear however 
that the infinity we need to consider is different from the one we were discussing to establish the 
logarithmic scaling of the energy. 

Near-flat space regime 

The near-flat space regime shows up when considering the large momentum limit of the dispersion 
relation (|5.37|) . Namely, keeping the leading contribution at p ~ cxd order by order in the 1/g 
expansion ()5.37p . we find thaO 



1 p"^ 7P'' 
2p Ag g^ 



E-P= — -TZ + \ + --- ^ (5-41) 



with 7 a numerical constant, see Eq. (|5.38|) . We immediately recognize the pattern observed previously 
for scalar. Namely, the suppression by a power oi 1/g comes with an enhancement by p"^. The near- 
fiat-space regime resums these contributions at fixed energy E = E/g^/^ and momentum p = pjg^l'^ . 
In this scaling limit, the dispersion relation ()5.4ip reads 

E-p = bEl^g^O{\lg), (5.42) 

with the correction bE determined at small p by ()5.4ip . In terms of the rapidity x, the near-flat space 
regime p ~ g^l'^ corresponds to working around ~ 1. More precisely, using (j5.39p . we get 

X- g = ^bx{p)^ ... , (5.43) 

where we have rescaled to the original rapidity, x = gx, and introduced 5x{p) = —l/2p^ + 0{p^). For 
the Bethe rapidity, u = x + g'^/x, it means that u — 2g = 5x{p)'^ + .... As we vary p from to -|- 00, we 



*°One needs also to include the subleading classical contribution, yl+p^ = p + l/{2p) + . . ., which turns out to be 
of the same order as the loop corrections in the near fiat-space regime. 
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expect that 6x{p) will range from — cxd to +00 and in particular changes sign. At some intermediate 
value of p the quantity dx{p) would vanish meaning that we are at the boarder between small and 
large fermion domain. It corresponds to the case where u — 2g is minimal and equals to zero. 

If we come back to the general representation (|4.34|) . the near-flat space regime u = 2g + 0(1) 
corresponds to expanding the function 7^(t) at g ^ 1 with t/g kept fixed. It is known how to 
construct the solution to the BES equation in this regime [64t [55] and an expression for the correction 
5E can be derived in principle for any p, at least parametrically. It would be interesting to perform 
this analysis in detail and check the matching with the giant hole regime when p^ \. 

5.3 Gauge field and bound states 

For gauge field and bound states, the string perturbative regime, E,p = 0{1)^ corresponds to taking 
the strong coupling limit while keeping the rescaled Bethe rapidity u = u/{2g) fixed and such that 
V? < 1. We will study this regime up to 0{l/g'^), but before going into detail it is interesting to see 
why this regime requires < 1. To this end, let us rewrite the relation ()4.44p . giving energy and 
momentum of gauge field {i = 1) and bound states (i > 1), as 



where e = i/ (4(7). When the coupling constant gets large, we have e — )■ and the contour of integration 
in (|5.44p becomes small, leading to 



in a first approximation. So the smallness of e has compensated the large factor of g, and energy and 
momentum are of the order of the function A(2gu), which is of order 0(1) as we shall see. We stress 
that this is correct for < 1. Otherwise, when > 1, since the integrand A{2gv) in (|5.44p has 
cuts running from v = ±1 to ±00, respectively, the contour of integration does not become small: it 
wraps around one of the cuts, and energy and momentum are of order 0{g). This is the giant hole 
regime |15] . whose analysis will not be given below because to leading order the dispersion relation is 
the same as for scalar and fermion, as said before. The transition through the points li = ±1 requires 
a resummation of the perturbative expansion that can be done by working with u^2g = 0(1). This 
is the near-flat space regime that we will consider mostly to leading order. Finally, we will consider 
another semiclassical limit reached when £ ~ (7 ^ 1. In this limit, we have e = 0(1), meaning that the 
contour of integration in (I5.44|) does not shrink at strong coupling, and both energy and momentum 
are then of order 0{g). Physically, this limit amounts to giving to the GKP string a (macroscopic) 
transverse spin S± = i ^ g. We will derive the leading order expression for the energy at arbitrary 
e, assuming an excitation at rest, p = 0. 

Perturbative regime 

As said before, in the perturbative regime, g ^ 1, u = u/{2g) fixed and u'^ < 1, energy and momentum 
are 0(1). Starting from ()5.44p and performing a bit of algebra, as explained in Appendix (DJ one 
observes that the dispersion relation can be written as 




(5.44) 



E + p 



YA(2<7u) , 



(5.45) 




(5.46) 
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where /(u) admits a 1/(7 expansion at fixed u From the knowledge of f{u), one can easily recon- 
struct both energy and momentum by separating odd and even parts in u, 



E 



i 


' n + u 


7^ 


\l-u 


i 


'n + u 


72 


\l-u 



) fin) + 

'4 



1 



1/4 



1 + U 

l-n\i/4 



1 + u 



fi-u) 



(5.47) 



For the dispersion relation, it is more convenient to observe that 

2ff{u)f{-u) . 



(5.48) 



For instance, using the leading order solution to the BES equation [62j . one immediately finds that 
(see Appendix dD]) for details) 

/(n) = 1 + 0(1/5) , (5.49) 



giving 



E^ 



P 



2t 



(5.50) 



This is a relativistic dispersion law for a particle with mass \f2i. In particular, it is saying that the 
gauge field {i = 1) has mass \/2 at strong coupling. There are exactly two bosons with precisely this 
mass in the spectrum of string fluctuations [13j, that we verify to correspond to the two twist-one 
excitations Fj^^\±i2 of the gauge theory. The linearity in I of the mass formula indicates that a bound 
state is no more than a pile of i gauge fields at strong coupling. This is presumably a consequence 
of the fact that there is no interaction between them at this order. As we shall see the masses still 
satisfy a linear law at one loop and the gauge fields get bounded starting at 0(l/g^). 

Including radiative corrections to the dispersion relation (j5.5Up can be done order by order in the Xj g 
expansion with the help of the results of [Ml EH ES] for the solution to the BES equation. The details 
of this analysis for corrections up to 0(1/^^) are given in Appendix [Dj The dispersion relation so 
obtained is simplified when expressed in terms of rescaled energy and momentum, Er = E/i,pr = p/i- 
For further simplification, we will simply omit to make use of the subscript r in the following, keeping 
in mind that the correct expressions require the substitutions E ^ E/i,p ^ p/i. After these brief 
remarks, using the expression quoted in (jP.lSp for the the function f{u), the dispersion relation (I5.48P 
is found to be 



E 



2 + p2 



{l+p^y (a + (3p^ + jp^ + 5p^){l + p 



2\2 



85(2 +p2; 



+ 



ff2(2+p2)2 



+ Oil/g' 



(5.51) 



with 
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The {u,p) relation reads, up to 0{l/g), 



u 



pV2 + p2 
1 -Fp2 



1 + 



(3 = - 
5 = - 

a + hp^ 
5(2+ p2) 



'48 ' 
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2^ 



(5.52) 



+ 0{l/g') 



with 



a 



1 _ 3 log 2 
8 2^ 



1 _ 3 log 2 
8 4^ 



(5.53) 



(5.54) 



*^Note that, in spite of the notation, the function f{u) in (|5.46p is different from the one entering the relation H5.32p 
for the fermion. Their strong couphng expansions are however very similar when both expressed in terms of u. 
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As a particular application of <\5.51\ [532|l . we find the (rescaled) mass 



m 



E{p = 0) = V2 



1 

16g 



a 



(5.55) 



We note that the first correction to the mass is negative, indicating that it decreases with the coupling. 
This is in qualitative agreement with the expectation that the mass increases monotonically between 
its weak and strong coupling values, 1 and ^J2, respectively. Next, recalling that the mass m in (|5.55|) 
stands for a rescaled one, we observe that the breakdown of the linearity in I appears at Oil/g^ ) and is 
controlled by the term ~ — in a, see Eq. (|5.52|) . Due to the negative sign of the latter contribution, 
we conclude that the gauge fields are effectively bounded by the interaction. It is curious nevertheless 
that the latter phenomenon is delayed to 0{l/g'^), and it would be interesting to look at the scattering 
between gauge fields at strong coupling to shed light on it. 

Examining futher (j5.5ip we find singularities at = —2. They are of the same sort as the ones 
encountered for scalar in the non-perturbative regime. They account for renormalization of the mass 
and one can get ride of them easily by rewriting ()5.5ip as 



E = \l m? + 



1-^ + 



p2(/3' + 7y + ,5/) 



9' 



+ 0{l/g^) 



(5.56) 



with m given in (|5.55|) and the coefficients /3','y', found as 

3 log 2 
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(5.57) 



The expression (|5.56|) is obviously simpler than the previous one (|5.5ip . However, it obscures that 
the dispersion relation satisfies the identity |14j 

(5.58) 



Eip = ±i) = l, 



which is apparent in (j5.5ip . This identity can be proved to be exact. The proof will not be given 
here and we simply note, using (|5.53p . that the points p = ±i corresponds to u = zkioo. Notice also 
that after rescaling to original units, the relation (I5.58P reads E{p = -^itj = i, with E the energy of 
a bound state of £ gauge fields. 



Near-fiat space regime 

The strong coupling expansion of the function f{u) that gives energy and momentum ()5.46p is singular 
around the points v? = 1, see Appendix |Dj Around u = 1, say, it turns into an expansion in inverse 
powers of 2g(l — u) = 2g — u indicating that a resummation is needed. Using the expression ()5.53p . 
we note that the point u = 1 requires taking the momentum p to be large. More precisely, the 
asymptotics at large p reads 

^ = 1-^ + 0(1//), (5.59) 

to leading order at strong coupling. Thus working around the points u ~ 2(7 is equivalent to consid- 
ering momentum p ~ g'^^^- This is the near-fiat space regime. From the result (j5.5ip we find that 
the energy E = Ejg^l"^ is given by 

^ = p + 5^(j5)/^ + 0(1/5), (5.60) 
where p = p/g^^^ is kept fixed and with 5E{p) determined at small p by (|5.5ip as 

SE{p) = \--p^ + 6p'^ + 0{p^^), (5.61) 
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with 5 given in ()5.52p . We note that to leading order E = j5+0(l/y^), even after the rescahng to the 
natural units {E,p) — )■ {E/£,p/i). We verify thus the flavor independence of the leading dispersion 
relation in the near-flat space limit. The subleading correction 5E{p) does depend on i however, via 
the coefficient 5 in (|5.61|) . see Eq. (|5.52p . 

Note finally that it is possible to derive a parametric representation for the subleading dispersion 
relation valid for any p. It is however difficult to deal with it and the matching with the giant hole 
domain is not obviously observed. It would be interesting to work it out in detail. 



Semiclassical bound state 



As mentioned before, it is possible for bound states to reach another semiclassical regime without 
taking the momentum to be large. We can assume p = and choose i ^ g- It corresponds to giving 
a macroscopic transverse spin S_i_ = £ ^ g to the GKP string. For p = the energy is minimal at 
given spin S±. This regime is easy to study to leading order at strong coupling, the mass of a bound 
state reading 

iv J 



E{0) = 2V2g J 



dv 



1 



(5.62) 



with e = l/Ag = S±/4:g. The relation (|5.62|) is a direct consequence of the general formula (|5.44|) 
given the strong coupling expression for A{2gv), see Appendix[Dj Note also that energy in (|5.62|) has 
not been rescaled by i, in distinction with the convention of the previous subsection. The integral 
in ()5.62p can be taken exactly yielding 



E{0) = 4g ( Ji]{r, - 1) + log r]^/^ + ^/^f^ 



0{g' 



(5.63) 



where rj = vT+~i^. The mass is thus a smooth, monotonic increasing function of e, interpolating 
between 



at e ~ 0, and 



E{Q) = AV2g 



e3 



+ 0{g^) = V2Si_ + ..., 



(5.64) 



E{0) = 2g 



2e + (log (4e) - 1) + 0(l/e) 



+ 0(^0) = 5x + 2g (log 



Si 



1 + 



(5.65) 



at e ~ cxD. It is interesting to note that the quantity E{0)/{\/^S±) is always smaller than 1 and 
decreases when S± increases. We verify again that the mass of the bound state is smaller than the 
sum of the masses of its gauge-field constituents. 

Soliton above the GKP vacuum carrying spin S± ~ g can be found directly in string theory, and, 
remarkably enough, has energy at rest exactly given by dSSSDlll This soliton should also be some 
special case of the two AdSs-spins solutions constructed in [67]. In the limit S± = £ ^ o, i.e. e ^> 1, 
the asymptotics behavior (j5.65p is also found to be in agreement with the result of 



6 Conclusion and outlook 

In this paper we have classified excitations propagating on top of the long GKP string and studied 
their dispersion relations both at weak and strong coupling. Developing ideas of seminal analy- 
ses [351 El EEl El dH] we have derived large-spin integral equations from the all-loop asymptotic Bethe 
ansatz equations for the spectrum of scaling dimensions [2J. The equations devised in the current 

am grateful to J. Maldacena for sharing with me his notes on this solution. 
*^In the notations of [2Uj . the limit e ^ 1 should correspond to a — >■ cxa (and L — 3), with a — S/£ and with S the 
total spin. Note also that our result gives the GKP energy of the semiclassical bound state measured above the vacuum. 
So one should substract the vacuum contribution in the result of [50] to make the comparison. 
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consideration are of the type encountered previously in the Uterature [IHl EJ [l9]. There exists an 
extensive methodology for their solution [STJ [USl ESI [Ml ESI ES]- In particular, owing to their close 
relationship with the BES equation, which controls the vacuum distribution of roots, we were able to 
construct an exact parametric representation for the dispersion relations, which, though not totally 
explicit, only requires the knowledge of the solution to the BES equation. This reduction in complex- 
ity allowed us to analyze the dispersion relations in various regimes both at weak and strong coupling. 
As particular applications, we have derived the first few terms in their weak coupling expansion and 
computed corrections to the spectrum of string fluctuations at strong coupling. A natural next step 
is a thorough analysis of the complex momentum plane for all excitations, and, in particular, of the 
Wick rotation symmetry proposed in [21j . 

At a more dynamical level, we have identified the mechanism that leads to the restoration of the 
SU{2) symmetry accompanying the shift from the BMN to the GKP vacuum. We have also observed 
the enhancement of the residual symmetry group of the GKP vacuum in the weak coupling limit 
to /S'L(2|4). This peculiar dynamical aspect of fermions was related at the level of their dispersion 
relation to the necessity of working with two Bethe rapidity planes to cover their full kinematics. 

Among the dynamical issues that would be worth pursuing is the fate of the missing mass 2 boson. 
There are possible decay channels for it, in two fermions, for instance, but it would be enlightening to 
give a more quantitative estimate of their amplitudes. One might want to understand, for instance, 
whether the mass 2 boson is stable to any order at strong coupling and gets destabilized by non- 
perturbative correctionslll The question of stability of the mass 2 boson is also related to the 
completeness of the basis of excitations that have been identified so far. Another dynamical issue 
concerns the size of the wrapping effects. They were assumed to be associated with 0(1/5*) corrections 
based on physical expectations and direct computations in the s[(2) subsector [52]. The latter result 
only covers the case of scalar excitations and it would be interesting to check whether the estimate 
holds true for other flavors as well, by means of the Y-system proposed in [^. 

There are possible applications of the analysis performed in this paper to some related problems. 
One of the most interesting one concerns the computation of scattering amplitudes along the lines 
of [21j . Another possible application is to the construction of the spectrum of excitations around the 
GKP string embedded into the string dual of the Aharony-Bergman- Jafferis-Maldacena theory ^68j . 
The all-loop asymptotic Bethe ansatz equations are known for this theory as well [691 [70] and are 
closely related to the ones for the planar = 4 SYM theory. From their structure we expect that the 
spectrum of large spin excitations can be analyzed along the lines presented here. Moreover, up to 
the replacement of the coupling constant by the interpolating function [71^ i70j of the Gromov-Vieira 
equations, the dispersion relations should be directly expressible in terms of the ones constructed 
here. 
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A Basic integrals 

In this appendix, we provide basic integrals, mostly borrowed from [18| . that are useful to performing 
the all-loop analysis given in the Section O 

''''This seems a quite reasonnable expectation since its presence in the spectrum of string fluctuations is tied to the 
breakdown of the SO{6) symmetry down to SO{5), which is only restored at the non-perturbative level. 
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Basic one-loop integrals 

Expressions of the type 

^ , , , fr + 2iu\ f2u\ , , , 

1t{u) = — zlog = zarctan — , (A.l) 

\r — 2m/ \ T J 

with asymptotics ibvr at u = ibcxD, respectively, make their appearance below, after separating the 
one-loop scattering phases from their radiative corrections. A convenient representation can be found 
as 

POO rjf 

I,{u) = 2 — sin {ut) e-''*/^ . (A.2) 
JO t 

Basic all-loop integrals 

The building blocks for the formulae given below are the following integral representations for inverse 
powers of the shifted all-loop spectral parameter x. They read 



(TiT ^J„(2<7t)e±-*-P*/2, (A.3) 
Jo t 



where x^^^^u) = x{u ± ip/2) and x{u) = {u + \Ju^ — {2gY)/2. Here, the index n belongs to N* and 
the function Jn{t) is the n*^ Bessel's function. When p = 0, one can also find the representations (for 
> {2gf) 



^\ = (-) = (-l)"(2n) / - J2„(2gt) cos {nt) , 
xj \gj Jo t 



(-l)"-i(2n - 1) / -J2n^i{2gt) sin {ut) 
Jo t 



where x' is the second solution to u = x + g'^ /x, reading x'{u) = (n — \J — {2gY)/2. The identi- 
ties (1X41) are consequences of ()A.3P when taking into account that 

/ -J2„(2(7t)sin(ut) = 0, / — J2n-i{2gt) cos {ut) = {) . (A.5) 
Jo t Jo t 

for > {2g)'^ and n > 1. 

Spin-chain higher conserved charges 

With the help of ()A.3p . the spin-chain higher conserved charges 



where x^^ = {x^{u))^, can be found as 



2( — \] f°° at 2( — 1) f°° at 

Q2n = ^ J , / -J2n^l{2gt) COS {ut) e-*/^ , Q2n+1 = -^-^ / —J2n{2gt) siu {ut) e-*/^ 

5^""^ Jo t g^^ Jo t 

(A.7) 



with n > 1. 



Integrals for inhomogeneous terms 



The all- loop large spin equations considered in Section [3] contain various inhomogeneous terms. Useful 
integral representations for them can be obtained with the help of few equalities. 



50 



The first identity holds for 



, / x-\ , fl + 2iu\ , [l + g'^/x+^\ , 
j(x) ^ _,iog -— = -zlog -— - - zlog , , , (A.8) 



x+ j \1 — 2iu J \l + g'^/x 

that involves the shifted spectral parameter = x{u it i/2), and where = (x^)^. By Taylor 
expanding the right-hand side of (|A.9P around g = 0, and making use of the integral representa- 
tions ()A.2l IA.3P , one obtains 

I{x) = 2 — Jo(25t)sin(ut)e"*/2 _ (^ g) 

JO t 

The second equality holds for 

1 - g'^/x^y^ 
1 - c/Vx-y[+Pl / ' 



/p(x, y) ^ -i log ( , ) , (A.IO) 



that involves the two shifted spectral parameters x^ = x{u ± i/2) and yt^Pl = y[v it ip/2). Going 
along the same lines before, one gets 



(A.ll) 



roo rOO 

Ip(x, y) = {2gf / dt cos {ut) e-*/^ / dsK{2gt, 2gs) sin {vs) e'P''/'^ 
Jo Jo 

poo roo 

- {2g)^ / dt sin {ut) e'^^^ / dsK{2gt, 2gs) cos {vs) e'P'/^ , 
Jo Jo 

where the symmetric kernel K{t, s) = K{s, t) is given in (j3.23p . In the particular case p = 0, we have 



Io{x,y) = -ilog = 77— , (A.12) 

\ 1 - g^/x y I \1- y'/x I 



with y = {v + \/v'^ — {2g)'^)/2 and y' = {v — y/v"^ — {2g)'^)/2. It can be evaluated as 

/■OO poo 

Io{x,y) = {2g)'^ dt cos {ut) e'^/"^ dsK+{2gt,2gs) sin {vs) 
Jo Jo 

poo poo 

- {2g)^ dt sin {ut)e-^/^ dsK^{2gt,2gs) cos {l 
Jo Jo 



(A.13) 



(VS) . 



for v"^ > {2g)'^ and with the kernels K±{t, s) given in (j3.24p . The last two relations are easily deduced 
from (lAlTl [3:231 irap . by applying fO) . 

Finally, to obtain the inhomogeneous terms for the gauge field and its bound states, we need to 
consider 

We find that 

-ps/2 



roo roo 

fp{x,y) =2{2g)'^ dt cos {ut) e'^/"^ dsK+{2gt,2gs) sin {vs) e 
Jo Jo 

poo poo 

-2{2g)^ dt sin {ut)e~^/^ dsK^{2gt,2gs) cos {vs) e 
Jo Jo 



(A.15) 



B Formulae for all-loop analysis 

In this appendix, we explain how to derive the expressions for ^dressing (^); -^dressing (^); Starting from 
the representation given in ^ for the dressing phase. Then, for all of the momentum-carrying excita- 
tions considered in this paper, we construct the inhomogeneous contributions X*(m, Ui,),ii,{u, u*). We 
give also the defining relations for the auxiliary functions that permit to cast the integral equations 
into the form of infinite systems of equations. Finally, we provide expressions for the infinite series 
of source terms entering the latter set of equations. 
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Dressing phase 

According to [21 [551 [2], the dressing phase a{u,v) = exp {i9{u,v)) admits a bihnear expansion over 
the individual charges qr{u),qs{v), see Eq. ()A.6p . An exphcit expression for it was proposed in [3] 
and reads, after separating even/odd parts in u, 

roc 

e{u,v)\oddinu = 2y2{2n) / -J2n(25t)sin(nt)e-*/2 

^1 Jo s e* -1 

0(u,t;)|evenmn = -2 V(2n- 1) / — Jsn-l (2c/t) COS (ut) 6 



X V (-l)™2(2m),- / , 



5] (-l)™2(2m)<72™ / 
m>l ^0 



1 



with Jnit) the Bessel's functions. To get the expressions for Xdressing('u),^dressmg(^t)) we perfom the 
replacement v — >■ u^j in (jB.ip and then sum over all the main roots u^j. This has the effect 
of substituting qr{v) in (jB.ip by the total charges Qr, since the roots n4 are the only spin-chain 
momentum carriers. For our purposes, we only need to consider total charges measured above the 
vacuum, qr = Qr — (7™'^""™. Then, introducing the generating functions as in ()3.30p . it 

is not difficult to evaluate 

K4. 

(u) = 2du ^ 9{u, Uij) , (B.2) 

and to cast it into the form (|3.29p . 
Scalar 

By definition ()3.261[3r27p . the inhomogeneous terms Xh(u, Uh),2h('u, Uh), associated with a hole carry- 
ing a rapidity Uh, immediately follow from the knowledge of the undressed scattering phase between 
main roots. The latter quantity is given in ()3.1ip . Hence, combining ()3.26l [3?27l 13. lip and using the 
formulae ()A.2l lA^ lA.lip of Appendix [Al one easily finds that 



Xh(u, ?ih) = — 2 / dt cos {ut)e I 7rh(t, Uh)) 
2h(ii, Uh) = — 2 / sin (tit) e~*'^ 7fh(t, tih) J 



(B.3) 





with 

/■oo 

7rh(t, Uh) = cos (uht) e-*/2 -Joi2gt) - (2gf / ds tK{2gt, 2gs) cos (uhs) e'^/^ 

JO 

7fh(t, Uh) = sin (uht) e"*/^ -(2c/)^ / ds tK{2gt, 2gs) sin (uhs) e"''/^ , 

Jo 



(B.4) 



with the kernel density K(t, s) as in ([3231 [323) • 

The defining relations for 7^ (t) , 7^^ (t) , that permit to cast the integral equations (|3.34p into the 
form of the systems of equations (I3.40|) , are 

-tl{2gt) + -i'l{2gt) = (e* -l)Vt^{t) - cos {n,t) e-*/^ +jQ{2gt) , 
^l{2gt) + 7^(25t) = (e* -l)Oh(t) - sin (uht) e-*/^ . 
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The infinite series of source terms for the system ()3.40p read 



= -l ^'-^ (cos iu^t) eV2 -Joi2gt)) 
Jo t —1 



(B.6) 



Finally, evaluating (jB.SP at small t, we get the relations 

f^h(0) = -^ + 257?(5), Oh(0)=tZh + 257?(<7), (B.7) 
where the coefficients 7i,7i, originate from the small-t asymptotics 

7^(2fft) = 257it + 0{t') , 7-(25t) = 2gj\t + 0{t^) . (B.8) 
Note that we used the fact that 'j\{t) ~ 7+(t) ~ t"^ at small t. 

Gauge field and bound states 

To compute the inhomogeneous terms (|3.26|) associated to a gauge field {£ = 2) or its bound states 
{£ > 1), carrying a rapidity ui, we need the expression for their scattering phases with a main root 
u. These excitations are embedded into the asymptotic Bethe ansatz equations (|2.1|) in the form 
of stack of strings [20]. For instance, considering gauge field and bound states with positive so (2) 
charges, they are found as a stack of one {£ + l)-string of roots U3, one ^-string of roots U2 and one 
{i — l)-string of roots ui, with a (real) center-of-mass rapidity ui. Then the scattering phase for a 
main root u going through the stack is obtained by taking the product 

£+1 t £-1 

Su{u,U£) = n -5*43(^^,^*3,^) n Si2{u,U2j) Y[ S4l{u,Uij), (B.9) 
j=l j=l j=l 

where 

U3,j-w='-{2j -2-£), U2,j - ue = '-{2j - 1 - £) , - = '-{2j - £) . (B.IO) 

The notations ^43 , ^42 , S41 , in (jB.OP stand for the elementary scattering phases encoded in the middle- 
node equations of the asymptotic Bethe ansatz equations (|2.1|) . They read explicitely as 

Si3{u,U3j) = ^, Si2iu,U2,j) = , 54i(m,U1j) = 5- ^. (B.ll) 

X- - X3j ' 1- g^/x-xij 

Plugging (|B.10p into (|B.9P and making use of (jB.lip . one finds after a bit of algebra 

a „ N _ ^ - + 2 1 - r/x x\ 'I- g'/x x\ 

^4£(^^, Ul) - [_^] , (B.12) 

u — ui 1 — g^ jx^Xm 1 — g^ lx'^x\ 

where x"^^ = x{u£ + i£/2). For later convenience, we have factored out the one-loop scattering phase 
in ()B.12p . which is the one expressed in terms of the difference of the Bethe rapidities u — U£. That 
was done by using algebraic identities between the deformed spectral parameters x, X£, and the Bethe 
rapidities u,ui (see [1]). The result ()B.12p is a straighforward extension to ui ^ of the findings 
of [2D]. 

Equipped with the expression (jB.12p . we can now evaluate the inhomogeneous terms Zi{u,Ui), 
ii{u,Ui), as defined in Using the formulae (1X21 [Xl |AT5]) . one can show that 



Ie{u,U£) = —2 / dt cos {ut)e ^^'^ Tri{t,ue) , 

^?oo _ (B.13) 

Xi{u,ui) = —2 / dtsin(ut)e vf^ (t, u^) , 
Jo 
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where 

roo 

7ri{t,ui) = cos (lift) e-^*/2 -M^gt) - {2gf / dstK^{2gt,2gs) cos {uis)e-^^/'^ 

Jo 

roo 

jf(,{t, Ui) = sin {uet) e'^/^ -(^2gf / ds tK+{2gt, 2gs) sin (ues) e'^^/^ _ 

Jo 

with the even/odd kernel densities K±{t,s) as in (j3.24p . 
The defining relations for 7j_(t), 74(i), are 

^i{2gt) + Y_{2gt) = (e* - cos (ut) e'^/^ +Jo(2(7t) , 

For the two infinite series of source terms in (|3.4U|) . we found 

4niue) = -f ^l^f^ (cos {u,t) e-«/2 -U2gt)) , 



(B.14) 



(B.15) 



(B.16) 



and 



-£ . ^ r dt J2n-li2gt) . _ft/2 

i^2n- iiui) = - — — j —T- sm e ^ 

Jo t 1 — e ^ 



(B.17) 



Finally, we note that 



nei0) = -- + 2gji{g), ne{0) = ui + 2gji{g), (B.18) 



where the coefficients 7i,7i, originate from the small-t asymptotics 

7i(25t) = 257ft + 0(t3) , ^^_{2gt) = 237ft + 0(t3) . (B.19) 

Large fermion 

For a large fermion, the physical rapidity X{ has modulus bigger than g and is given in terms of the 
Bethe rapidity nif as 

Xf=^(nif+^n2-(25)2) . (B.20) 

We use the subscript 'If (for 'large fermion') to keep track of this information when working with the 
Bethe rapidity. Note that for real rapidity Xf we need uf^ > {2g)'^. 

The phase for the scattering of a main root u with a fermion carrying a large rapidity Xf is conve- 
niently written as 

SaAu, Xf ) = f- . (B.21 

u - u\i- -^1 - g'^/x+Xi 

Given the expression ()B.2ip . we can evaluate the source terms X\i{u, tiif),Z'if(M, nif), as defined in (j3.26p . 
Using the formulae ([Ml Ell lAlUlAlH]) . we find 

roo 

T\f{u, u\i) = — 21 dt cos {ut) e~*/^ TTif (t, uif) , 

^'00 (B.22) 

iif (it, uif) = — 21 dt sin (lit) e vrif (t, uif) , 
Jo 
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where 



POO 

7rif(t, uii) = cos (nift) — Jo{2gt) — 2g^ / ds tK^{2gt, 2gs) cos (uus) , 

Jo 

POO 

Tf\{{t, uif) = sin (uifi) — 2(7^ / ds tK^{2gt, 2gs) sin (uifs) . 

JO 

with K±{t,s) given in ([231D- 
The defining relations for j^].{t) , j]^{t) , are 

jl{2gt) + 7-(2(7i) = (e* -l)nu{t) - cos (uift) - Jo{2gt) , 
jli2gt) + j^li2gt) = (e* -l)i7if(t) - sin (uift) . 

The source terms in (|3.40p are found to be 

If, , dt J2n{2gt) , , , ^, I dt ^ . , , , 

'«2n(«lf) = - / -f r- (^OS (^If*) ~ M^gt)) - o / —J2n{2gt) COS (Ulft) , 

Jo r e'^— 1 2 Jo t 

If / N dt J2n^l{2gt) J -^\\ 

Jo t e'^ —1 

and 

,lf f^dtJ2ni2gt) . 



(B.23) 



(B.24) 



(B.25) 



t e*-l 



«^2n-i(^lf) = - / -I -f — sm {ui(t) - - / —J2n-i{2gt) sm (U\it) . 

Jo t —1 2 Jo t 



(B.26) 



We note also that 



Oif(O) = 2g-ff , j^if(O) = uu + 2g^f , (B.27) 
where the coefficients 7^ , 7'^, come from 

j'li2gt) = 2gj{h + 0{t^) , j'l{2gt) = 2gj'^t + 0{t^) . (B.28) 

Small fermion 

For a small fermion, the physical rapidity Xf has modulus smaller than g and is given in terms of the 
Bethe rapidity Ugf by 

Xf = ^ (usi - ^ul-{2gY^ . (B.29) 

Same remarks as for large fermion (see below Eq. ()B.20p ). 

The phase for the scattering of a main root u with a fermion carrying a small rapidity Xf is 
conveniently written as 

54f(n,xf) = ^fl^^') . (B.30) 

Plugging this expression into (|3.26|) . we find that the inhomogeneous terms, which source corrections 
to the twist-two solution, are given by 

( \ — Xf/x+(n) \ 

Xsf(^i, ^isf) + 2sf(u, tisf) = -idu log , (B.31) 

\l- Xi/x (n) J 

where Xsf(n, nsf),2sf(^*, ^^sf), are even, odd, in u, respectively. Using the formulae (|A.121 IA.13p . the 
latter quantities can be found to admit the representation 

roo 

Isi{u, Usi) = - 2 dt COS (ut) e"*/^ TTgf (t, Usf) , 

7oo (B.32) 

Xsf (it, Usi) =-2 1 dt sin (ut) e vTsf (t, Usf) , 
Jo 
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(B.33) 



(B.34) 



where 

/■oo 

vTsf (t, Usf) = / ds tK^ {2gt, 2gs) cos (ugfs) , 

roo 

vfsf (t, Usf) = 2g'^ / ds {2gt, 2gs) sin (usfs) . 

The kernels s) are given in (I3.24p . 

The defining relations for ^^{t),^^{t), are especially simple. They read 

7f(25i)+7i'(25t) = (e*-l)l^sf(t), 
fl{2gt)+fl{2gt) = {e'-l)^,,{t). 

The non-zero source terms in (j3.40p are 

roo 

l<2niUsf) = - —J2n{2gt) COS (Usf*) , 

2 Jo t 

'^In-lC^sf) = o / —J2n-l{2gt) sin (Usff) . 

-ii JO I 

We have furthermore the relations 

Osf(O) = 2gfl , n,i{0) = 2gfl , (B.36) 
where the coefficients 7|^, 7^^, come from 

fl{2gt) = 257ft + 0(t3) , f_^{2gt) = 2g^ft + 0{t') . (B.37) 

Vacuum distribution 

The vacuum distribution density can be written in the large spin limit as 

/. ^ , /.sx ^ / , cos (ut) e*/2 Jo{2gt) - 1 



(B.35) 



2 2 /""^ 

Pvacuum('u) = - (log S - 1p{l)) + - dt ■ 

TT vr Jo 



vr 7o e* -1 (B.38) 

- (log5- -<5a0(n), 

where a^{u) = a^{—u) can be found as solution to the BES equation [3]. The function 5a^{u) = 
Sa^{—u) captures the subleading, 0{log^ S), correction to the vacuum (twist-two) distribution den- 
sity, considered in |53J. 

One can introduce the functions 0,^ (t) , 5il.'^ (t) , as the Fourier-like transforms of a^{u), 5a^{u), 
respectively, following (j3.3H [3T32|) . Then, defining 7^ (t), (57^(t), and 'y^ (t) , 5^^ (t) , respectively even 
and odd in t, by 

^f{2gt)+j^{2gt) = {e'-l)n<^{t), 
5jfi2gt)+5^^i2gt) = ie'-l)5n^{t), ^ ' ' 

one finds that 7±(t), ^7^ (t), solve the system of equations (j3.40p with source terms 



<5kS' = 2 



2g6n,i , 

- dt Jn{2gt)M2gt) - gt6n,i (^-40) 



t 



respectively. Note that our convention for introducing "y±{t) differs from the one used in [B^ by the 
overall factor ^4g\ Namely, we have 7±(t) = 4g'7-|-(t), where 7±(t) are the functions analyzed in [65J. 
The difference between the coefficients 6k^ given here and the ones used in [53] originates from 
the decomposition (|B.38p , which is different in the two cases. A redefinition of the type 6j^{t) 
57+(t) + a(l — Jo(t))i for some constant a, would lead to the coefficients of [53]- The reason for 
decomposing as in ()B.38p is that it improves the large n behavior of the coefficients 5k^. 
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C Isotopic SU{2) root 



In this appendix, we perform the all-loop analysis of the isotopic root u\, which implements the 
restoration of the SU{2) C S'[/(4) symmetry. We recall that this root is a stack formed out of one 
2-string of main roots, one root ^3 and one root u^, as 

n4 = Ub ± I , U3 = U5 = Uh, (C.l) 

where Ub is taken to be real|3 Our analysis also assumes that > {2g)'^ to avoid complication with 
the cuts structure of the all-loop asymptotic Bethe ansatz equations [2] . At the end the results show 
no dependence at all on this bound and can be easily analytically continued. It would be interesting 
nevertheless to understand how we should deal with the cuts structure to get explicitely the results 
when ul < {2gY . 

The appendix is organized as follows. We first prove that the stack (IC.lj) does not carry energy and 
momentum. The analysis also shows that they do not carry any of the higher spin-chain conserved 
charges To find the energy and momentum of the stack (jC.ip . we need to compute the correction 
to the distribution of main roots sourced by (jC.ip . This is done by contracting the scattering phase 
between main root and root Ub, from which an inhomogeneous term for the large spin integral 
equation is obtained. Then we show that the solution to this equation is exactly equal to its one-loop 
expression. This is enough to prove the statement that the stack is isotopic (jC.ip . Then, to establish 
that the stack (jC.l|) implements the restored SU (2) symmetry, we show that the effective Bethe 
ansatz equations for the roots Ubj are identical to the ones of an inhomogeneous SU{2) Heisenberg 
spin chain. The length of this spin chain is fixed by the total number of scalars (holes) and the 
inhomogeneities are given by the holes rapidities. 



C.l Energy and momentum 

The energy of the stack (|C.ip is the sum of the bare energy carried by the 2-string of main roots 
plus the energy induced by the deformation of the background. To get the latter contribution, 
we need the expression ph{u) for the correction to the distribution density sourced by (jClll . This 
function is given as the solution to the large spin integral equation. The main ingredient of this 
equation is the inhomogeneous term which follows from the undressed scattering phase between the 
probe main root u and the stack (IC.ip . The kernel of the equation is the same as the one already 
considered in Section O The main simplicity here, as compared to the momentum carriers studied 
in Section [31 is that the integral equation can be solved exactly for a stack. We recall finally that in 
our approach we consider separately the even and odd part in u of the correction density 

Ph{u) = -CFhiu) - ah{u) , (C.2) 

with the notations of Section [3J 

We start by constructing the inhomogeneous term. The phase for the scattering of a main root u 
through an isotopic root «b is given by 

S^hiu, -Ub) = 5*44 (u, u^'j Su (u, u^'j S43 (n, Ub) {u, u^) , (C.3) 

where 5*44, S43, 5*45, stand for the elementary scattering phases between a main root u and a root 
it4,it3,U5, respectively. Their expressions are given on the right-hand side of the middle- node equa- 
tions in (j2.ip . After isolating the contribution coming from the dressing phase, we find that the 
undressed scattering amplitude, denoted S^^, reads 



u 



Mb — ^^"""b + f / I — g'^/x^x- 




U-Ub + fu-Ub-^ \l- g-2/x-x[; ^ I 
*^This assumption can probably be relaxed to permit complex roots iibj- 

*®They contribute however to the total spin-chain momentum e'^ by a factor (—1)^'' where Kh is the total number 
of root Mb • 
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where = x{ua ± i). The phase ()C.4p generates the inhomogeneous terms Xb(u, nb),2b(n, Ub), 

that are given by 



i5„log54b(u,Ub) =Xb(n,'Ub) +Xb(n,nb) 



(C.5) 



where Xb(M, nb),Xb(M, Ub), are even, odd, in u respectively. Evaluating X\,{u,Uh),'ih{u,Uh), with the 
help of formulae given in Appendix [XJ we find 



2b(w,iib) 



dt cos {ut)e ^/'^ ■Ka{t,U\j) , 

) 

dt sin {ut) e^*/^ vf(j(t, Ub) , 



(C.6) 



where 



7rb(t,Mb) 

7fb(*,'"b) 



(l-e cos (ubt) + (2g)^ y dstK{2gt,2gs) cos {ui,s)e~ 
(l - e"*) sin (ubt) + {2gf J ds tK{2gt, 2gs) sin {u\,s) e~ 



(C.7) 



The kernel K{t, s) is given in (I3.23p . 

Now the claim is that the solutions Jlb(i)) ^b(i), to the equations (|3.34|) sourced by the terms (|C.7p 
are 

r^b(i) = cos (ubt) e~* , ^h{t) = sin (ubt) e~* , (C.8) 

where we recall that Q-\^{t)^Cl\^{t), are the Fourier-like transforms (j3.31l [3TT5|) of a\^{u) ^ a\^{u) , respec- 
tively. The fact that the two functions (jC.SP solve the equations would be straightforward if it were 
not for the dressing-phase contributions in (j3.34p . So let us see if the two functions (jC.Sp annihilate 
these terms separately. For the dressing-phase contributions to vanish, we need all the conserved 
charges sourced by a stack to be zero. The latter charges can be found as sum of a bare and induced 
contributions. Evaluating the induced contribution associated to (jC.SP we find 



^induced 
H2n 



induced 
y2n+l 



2(-l)" 
„2n-l 



2(-ir 

r,2n 



dt 
T 

dt 

T 



J2„-i(2gt)na(t) 



J2n{2gt)na{t) 



'2n 



2n 



2n 



2n-l 



2n-l 



2n-l 



h2] 



[-2] 



(C.9) 



But these are exactly minus the charges for a 2-string of main roots. It means that induced and 
bare conserved charges cancel one another. The dressing-phase contributions vanish for (jC.SP and 
the functions (jC.SP are exact solutions. 

We have shown that a stack with rapidity Ub do not carry any of the spin-chain charges. Therefore it 
does contribute to the anomalous dimension and cannot carry GKP energy. Moreover the correction 
to the density sourced by Ub is one-loop exact, 



Ph{u) = -CThiu) - (Tb(^i) 



1 

vr Jo 



dt cos ((n - Uh)t) e . 



(C.IO) 



From the large u asymptotics of pb('u), we get the GKP momentum by looking at the coefficient in 
front of the leading ~ contribution, see Eq. (|4.3p . But, the density (|C.10P scales like ~ l/u^ at 
large u, meaning that the stack is exactly momentum-less, ph = 0. This concludes the proof that the 
stack (jC.ip is isotopic. 



C.2 Effective equations 

To get the effective equations for the stack rapidities u^j, we pick up a probe stack, with a rapidity 
tib,fc, that we transport once around the spin chain. Gathering the phases for the scattering with all 
the possible roots, composing a generic state, we get the equations by imposing the condition that 
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the total phase equals 1. The last identity correponds to imposing periodic boundary condition to 
the Bethe wave function. 

The total phase can be decomposed into a product of a bare and an induced scattering factor, 
giving the equations as 

^ ^ ^induced ^bare ^ (C.ll) 

The induced phase takes into account the scattering with the main roots, 

5ir'"^^'^ = n'5b4Kfc,U4,i). (C.12) 
i=i 

Note that we take product over all the main roots, real plus complex, thus including the ones that 
form part of the stacks libjj with j ^ A;0 Though it may seem more natural to define the induced 
scattering phase by taking product over the real main roots only, we found simpler to gather the full 
set of main roots into a single factor. The immediate consequence of this prescription is that the bare 
scattering of the stack Ub j with its congeners boils down to a product of scattering phases with a set 
of fermionic roots carrying rapidities usj = = "Ubj. Since there is no bare interaction between 
fermionic roots, the only phases come from the scattering of the 2-string u^^ = -Ub.fc ± i/^ with the 
roots Us J = = tib,j, with j ^ k. Using p.ip we get 

s^r = \[ — ^1 , (C.13) 

r_|_Q] 

with a^b ^ = 2;(tib,fc i i). The equations are thus given by (|C.lip with the induced scattering 
phase ()C.12p and the bare scattering phase 

where the first term on the right-hand of (|C.14p comes from the propagation of the 2-string around 
the spin chain with length 2 + M@ while the last term is for the bare scattering with a generic set of 
M momentum-carrying excitations. With the help of the expressions computed below, we find that 
the equations (|C.lip read 

-rj Uh,k - Uh,j + 2 ^ TT Uh,k - Uh,j + I 

They are identical to the equations of a su(2) inhomogeneous spin chain of lenth K\i, with inhomo- 
geneities given by the set of holes rapidities u^j- We recall that K^^ is the total number of scalars, 
whatever su(4) polarizations they carry. Including scattering with the su(2) x su(2) C su(4) isotopic 
roots U2,ue, the equations (IC.ISP become the middle-node equations for an inhomogeneous 50 (6) spin 
chain, as follows from multiplying the right-hand side of ()C.15P by the factor 

K' ■ K' 

-Q ^b,fc - U2,j - § j| Mb.fc - uej - I 

j^l Uh,k — U2,j + § Ut,,k — U6,j + § 

where K2,Kq, are the numbers of solitary roots U2,uq (i.e., those roots which are not part of some 
stacks). The simplicity of the expression ()C.16P comes from the absence of higher- loops corrections 



*^We should not include scattering of a stack with itself, hence the condition that any product over the set of stacks 
should be restricted to j ^ k. 



48 



See Eq. (|2.16|) and comments following it for the origin of this quantity. 
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to the second and sixth node asymptotic Bethe ansatz equations ()2.ip . combined with the fact that 
there is no induced scattering phase since there is no couphng between roots of type 2, 6, and main 
roots. The simpHcity of (|C.15p . however, is due to a wealth of cancellations between the induced and 
bare scattering phase. To prove (jC.lSp . we will now compute S\^f^'^ and S™*^"'^'''^, starting with the 
former. 

Bare scattering 

We list below the various bare sattering phases 5'{jJ'^° for -k = h.,{,{,£,£. A simplification comes from 
the observation that the isotopic stack is invariant under the replacement {ui,U2,U3) — ^ {uj,ug,U5) 
changing the sign of the so (2) spin. Therefore, it is not necessary to distinguish between f, i, and f, i. 

• There is no bare scattering with a hole, 

Str = l- (C.17) 

• The bare scattering of the stack with a fermion reads 

[+2] _ 

5br = ^|^^V^. (C.18) 

• The bare scattering of the stack with a gauge field {£ = 1) or bound states {i > 1) is given by 

n^,k-u^-^{l+£/2) 1^ 1 - 5 V^g' 4'''' 1 " 5 V^ff 4"^ ' 

It is obtained by taking product of elementary scattering phases between the isotopic stack and 
the £-stack. 

This is the complete list of bare scattering phases for excitations considered in this paper. 
Induced scattering 

The induced scattering phase ()C.12p is more difficult to estimate, since it requires integrating over 
the continuum of real main roots. The starting point is the phase for the scattering of the root tib,fc 
with a main root u, 

ShA{uh,k, u) = S^4{uh,k,u)S^f'''''''^{uh,k, u) , (C.20) 
where for convenience we have factored out the the dressing phase contribution that reads 

^dressing ^^^^^^^^ = a"^ {u^,^, u)a'^ {u^,,, u) . (C.21) 

The undressed scattering phase S^^{ui^^k,u) is given by S^^{ui^^k,u) = S2^^{u,ui,^k), which, after 
using (|C.4p . yields 



CO/ ^ Uh,k-u-^Ub,k-u+^ (l- g /x[Jx \ 

5'b4(ub,fc,u = — ^— f . C.22) 

Uh,k-U+^Uk,k-U-^yi-g2/x\;]^X+j 

To proceed further it is useful to take the log of (|C.20|1 . We find that the latter quantity can be 
written as 

roc ij^-f- roo 

- ilogShi{u\,^k,u) = — sin (ub,fct) e"* u) - / — cos {uh,kt) 4>{t,u) , (C.23) 
Jo t Jo t 
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where 

<P{t, u) = 4>°{t, u) + u) , 4>{t, u) = 4>°{t, u) + 4>''~^{t, u) , (C.24) 

reflecting the decomposition (|C.20p . After a bit of algebra, using identities given in Appendix |Al we 
get 

roo 

cj)°{t,u)= - 4sinh(t/2) COS (lit) - 2(2^)2 / dstK{2gt,2gs) cos (us) e , 

Jo 



poo 

if[t,u) = - Asmh{t/2)sm.{ut) -2{2gf / dstK{2gt,2gs) sin {us) e 

JO 



■s/2 



(C.25) 



and 

/dressing/^ „ n ^ 9r9„\2 



(t,n) = 2(2g)2 r dstK^{2gt,2gsf-^^^^ 
Jo e^-1 

= 2(25)2 r dstK+{2gt,2gs)^±^^^ 
Jo 



(C.26) 



The kernels K{t, s), K^{t, s), K-{t, s), are given in (j3.23l [HT^^ . In analogy with the analysis in Sec- 
tional we have introduced generating functions q-{t, u),q^{t, u), for even, odd, charges q2niu) , q2n-i{u) , 
respectively, as 

= 2 ^(2n - 1)(-1) V"~'92n(n) Jsn-iW , 

(C.27) 

q+it, u) = 2 ^(2n)(-l)"5'>2n+i(n) J2n(t) • 



n>l 



The expressions above hold for a generic main root u. To obtain the complete induced scattering 
phase, we substitute u — )• u^j in (1C.23[) . sum over the set of main roots and take the large spin limit. 
Performing these steps leads to the decomposition 

J2 <P°{t, u,j) = J dv p{v)cp°{t, v)-Y, 0°(t, Uhj) + ^ cP°{t, u^j + i) + E '^°(*' «b,i - i) , (C.28) 

i=l j=l j^k j^k 

and similarly for Y^^=i 0°(^)^4,i)- Note that we took care of isolating the contribution from the set 
of 2-strings of main roots, Ubj i since they are not part of the continuum of real main roots. 
We also substracted the interaction of the stack Ub^fc with itself by restricting the sums to j ^ k in 
the right-hand side of ()C.28p . We could perform the same steps for the dressing phase contributions, 
but it is more convenient to note that since ^dressmg^ ^dressing ^ g^^.^ gjyen in terms of the generating 
functions for the individual charges, summing over the main roots amounts to replacing the latters 
by the generating functions for the total charges. Then we can write 

Ki M 

q.{t, u^j) = g— (t) + q_{t) , q.{t) = ^ q*_{t) , (C.29) 
j=l *=i 

where q-{t) is the generating function of the even charges measured above the vacuum. A similar 
identity holds for the generating function of odd charges, with q^^^^^(t) = 0, and we recall that 
there is no contribution to q- (t) , q+ (t) , from a stack u^j ■ 

It remains to take the integral over the density p(v) in (IC.28j) . In principle we should consider 
separately the contribution from the vacuum and the ones from the excitations. But there is actually 
a trick to unit the two. It turns out to be possible to trade the vacuum for a pair of two holes carrying 
the rapidities nh,± = iS*]!^ More precisely, the identity 

Pvacuum 

(n) = ph{u; Uh,+) + Phiu; tih,-) (C.30) 



''^This is reminiscent of the two large holes, with rapidities ~ ±5'/\/2, that appear in the Baxter approach [T7j. We 
stress however that here we need to choose the rapidities as ith,± = ±5* to get the right result up to subleading ~ log" S 
contributions. 
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is valid for both the ~ log S and ~ log*^ 5 parts of the vacuum distribution density, where ph{u', ^h,±) 
is the correction to the density sourced by a hole carrying rapidity Mh,±If^ We will not prove this 
relation here but we note that it generalizes the observation that the energy of a hole with a large 
rapidity Uh is exactly half the vacuum energy if Uh = ±5, see Eq. (|4.5U|) . Accordingly, to simplify the 
analysis, we will set the vacuum density to zero and work with an extended set of excitations. The 
two extra holes will be given the large rapidities tth,± = ±5* at the end. This trick can be apply to 
the charges in (|C.30P as well, such that we can completely ignore the vacuum contribution. 

We come to the most remarkable simplification of the computation of the induced scattering. 
Namely, we notice that after integrating with the density p{v) = —(t{v) — &{v) in (fCl28]) and adding 
the dressing phase contribution, we get precisely the left-hand side of the large spin integral equation, 
see Eqs. (|3. 341 [HTM]) . Using this equation then leads to 

K4, M+2 i^h ii'b 

j=l *=1 j=l j=l jj^k j^k 

(C.31) 

and similarly for J2f=i4'{'tjU4,j). The inhomogeneous terms 7r^(t, n,^), 7rb(t, 'Ub), are given in Ap- 
pendix [B] for the momentum-carrying excitations and in (|C.7|) for the isotopic stack. It is now 
straightforward to compute the induced scattering phase since all terms in ()C.3ip are known ex- 
plicitely. The result is 

Kb ■/l_n2/^[+2] \^A/+2 

^induced ^ TT '"-b.fc " Uh,j " ^ ^ /-^b,fc -^b,j -pj- ^induced 32) 

where 

• for a hole (scalar) with rapidity Uh, 

Cinduccd _ / ^h,k \ Uh,k — Uh .j, fn QQ\ 

Uh,k -Uh + 

for a fermion with rapidity Xf, either in the small or large rapidity domain, 

h2] \ [-2] _ 

cinduced _ | •"b,k \ ^h,k „ ,x 

•^bf - I -IZil [+2] _ ' (,<^.cS4j 

,k I ^h,k ~ 

• for a gauge field {£ = 1), or bound states {i > 1), with rapidity u^, 

' " l-K' ) - «« + 'd + V2) [ 1 - sv-t^r" 1 - 9V<^r" I ^ 

We observe that a large part of the induced scattering phase (|C.32p is cancelled by the bare scattering 
one. After simplification, we find 

A'h _ _ i Kb _ ■ 

I — 5bare ^induced _ TT ^b,fc ^h,j 2 TT ^b,fc ^b,j + ^ 




^"^This claim may look suspicious given that we said before that an important difference between the vacuum density 
and the correction sourced by an excitation is that the former scales like ~ log (S/u) at large u while the latter 

is suppressed as ~ 1/u. There is no contradiction however but simply an order of limit phenomenon. Namely, the 
large u behavior p*{u) ~ 1/u is correct only if the rapidity it* of the excitation is kept fixed in the process (it is 
actually an asymptotics for u ^ However, if the excitation is sent to infinity ii* ~ S 3> 1 and if only the 

leading ~ log and subleading ~ log*^ contributions are retained, then the large u asymptotics is transmuted into 
~ log(u«/M) ~ log(5'/?i), in agreement with the scaling of the vacuum density. Note finally that we need to send 
a momentum-carrying excitation to infinity for this to apply and since the vacuum has twist two we need two large 
twist-one excitations, that can be choosen to be two holes with zero total momentum. 
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in agreement with (IC.15[) . Note that in ()C.36|) we have already sent to infinity the two vacuum- 
building holes, with rapidities Uh,± = ^S. 



D Strong coupling solution to the BES equation 

In this appendix we recall few results about the strong coupling expansion of the functions 7±(t) 
and T^{t), at fixed t. They are relevant to the computation of the dispersion relations for fermion, 
gauge field and bound states, in the (string theory) perturbative regime. We also provide expressions 
for the analysis of the dispersion relation for a scalar. All the results given below can be obtained 
from [63l EU [38l ES] , but we will mainly apply the representations used in |i65j . 



D.l Formulae for fermion 

The (even/odd) functions J±{t) can be found to admit the representatiorl^ 

7^(it) ^ if (it) + ij^iit) = h^{t)Io{t) + hf{t)h{t) . (D.l) 
The special functions /o,i(*) defined as 

vr j-i VI — ^/ 71" J-i Ki^ — v/ 1 + v 

and normalized as /n(0) = 1 + n with n = 0, 1. The factors /lo'i(^) given at strong coupling and 
fixed t as 

hf{-it) = it- a{g)t'/g - ic{g)r' / + 0{t^/g^) , 
h^^-^t) = -th{g)t/g + d{g)t'/g' + Oit^g') . 
The expansion coefficients read 



1 31og2 / K 91og^\ ^^,,,2. , . ^ 31og2 K ^^,,,2. 
13 K 31og2 91og22 ^, . , ^ K 31og2 ^, , , 

(D.4) 

where K is the Catalan's constant. 

For computing the fermion dispersion relation in the perturbative regime, one needs to evaluate 

E+p = l--\{u), \{u)= ^0(2gt)e»*=/ -7^(t)e»S (D.5) 
I Jo t Jo t 

with u = u/{2g) and u'^ > (2g)^ for a real-momentum fermion. The computation of (jD.SP is easily 
done, order by order in the 1/g expansion, by using (|D.ll ID. 2} ID.31 ID^ , and after noting the basic 
integrals 

poo ._ _|_ 1 \ 1/4 

ij dt/o(-ii)e*"* = 2-2 (^l^j , 



poo 

I / dt hi-it) e'^^ 
Jo 



1 + u \u — 1 



When analytically continued to complex u, outside the cut < 1, the expressions in the right- 
hand side of (ID.6P can be used to evaluate the dispersion relation, in the perturbative regime, for 
physical (rescaled) rapidity x = x/g = u — yju?' — \ with modulus smaller than 1. Note that the zero 



^^We recall that one needs to perform a rescaling to derive expressions used in this paper from the ones given in |65) . 
Namely, 7^(i) = ^g"i±(t\ where 7±(f) are the functions constructed in [65| . 
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momentum limit corresponds to x = and thus to tt = oo. For illustration, to leading order at strong 
coupling one gets that 

X{2gu) = 2-2[^^j =2-2(^^j +0{l/g), (D.7) 

leading to 

/I _|_ ^\ 1/2 

E + P=[y^) +0{l/g). (D.8) 

As E,p, are respectively even, odd, in x, or equivalently u, one can extract the dispersion relation, 
E = E{p), from representation like (jP.Sp . To include higher- loop corrections, one can take derivatives 
with respect to u on both sides of (jD.6|) in order to generate new identities, and then apply (ID.H 
ID. 21 ID.31 ID.4P to compute (jP.Sp . Doing so, one finds that the dispersion relation can be written as 

/I + ^ \ 1/2 

E + p=(---] fix), (D.9) 



1 — X 



where f{x) is given to two-loops by 



2x(ax + b(l-x)2) 4xHcxi2 + x + 2x^) + dil-xf{l-x + x^)) , , 3^ 
f{x) = l 7- ^272 + 271 ^2U + 0(1/5' , D-10 

with the coupling dependent coefficients a, . . . , d, as in ()D.4p . 
D.2 Formulae for gauge field and bound states 

To evaluate the dispersion relations for gauge field and bound states, in the perturbative regime, we 
need the Fourier transform of T^{2gt), 

Mu)= I ^r0(25t)e»*, (D.ll) 

for V? < 1, where u = u/{2g). It can be obtained from the analysis of |65j . For < 1, it is found in 
the form of a 1/g expansion taken at fixed u, and it reads 

A{2gu) = -^-^ (\±^y^\{u) , (D.12) 

with 

.(,) = l + ^ + J^^) +^^ + ^^ + 0(1/9'). (D.13) 

g{l-u) g{l + u) g^{l-uy g^[l + uy 

The coefficients a{g), . . ., admit & 1/g expansion starting as 

32 levr \ 6144 647r2 bUn^ ' ^i^>' 



(D.14) 



1 31og2 / 7 K 91og22\ 2^ 
^(^) = -35-^+ 1,-6144 - 64;^ -^12^; /9 + 0{l/g), 

35 5K 15 log 2 451og^2 . 
= -6144 + 64^^ - ^12^ + ^12^ + ^(^/^^ ' 

d{g)= -^ + ^-^-^ + 0(1/,). 
2048 647r2 5127r 5127r2 ^ 

Now, to evaluate the dispersion relation, we need to integrate the function K{2gu), 

E + p = mg dvA{2gv), (D.15) 

J u—ie 
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where e = £/(4g). At strong coupling, we have e — )■ and for < 1 the contour of integration 
in (jP.lSP becomes small. Hence, one can Taylor expand in e to get 



E + p 



y 



A{2gu) 



■k{2gu) + 0{l/g' 



96(7^ dv? 

Making use of ()D.12p one concludes that Eq. ()D.16P can be casted into the form 



E + p = V2£ 



1-u 



with 



/(^) = ^(^) - TT^TTT^-^ + 0(1/9') ■ 



(D.16) 



(D.17) 



(D.18) 



96c/2 4(1 - li2)2 

Applying the relations ()D.18l ID. 131 ID.14P to ()D.17p we find the dispersion relations given in (|5.5H 

EM}. 



D.3 Formulae for scalar 

The dispersion relation for a hole in the non-perturbative regime is given in (15. ip in terms of an 
infinite set of coefficients which depend only on the coupling. They admit the representation 



8\/2 



(4n + 1)tt 
(4n + 3)7r 



-(4n+l)7rg 



-(4n+3)7rg 



1 



1 



=Re 



°° dt e*(*-''/4) 



(4n + l)7r 4\/2 [Jo t + i{4n + l)'iTg 

°^ dt e*(*+''/'^) 



+ 



(4n + 3)7r 4^/2 



=Re 



t + i{4n + 3)7r5t 



(D.19) 



with T^{t) = r'^(— t). Note that despite their appearance, the relations ()D.19P do not assume that g 
is large, but only g ^ 0. At g = for instance, the integrals in ()D.19p are subleading, giving 



rrinig = 0) 



8V2 



(D.20) 



Here we are mainly interested in their expression at strong coupling. Going along the lines of 
we get the coefficients ()D.19P in this limit as 



m4„+i = (Sg) e-(^-+i)-^ r(f )(2(4n + Ijvr^)-^/^ 



r(n + i) 
r(n + i)r(i) 



1 + 



3 31og2 
32(4n + l)iT levr 



m4„+3 = (Sg) e-(^-+=^)-^ r(f)(2(4n + 3)vrg)-^/^ /^^" + ;^3 

r(n + i)r(|) 



/9 + 0{l/g') 
3 



(D.21) 



1 



5 31og2 
32(4n + 3)7r IGtt 



/9 + Oil/g') 



We recover for mi the expression of the mass gap of the 0(6) sigma model obtained in [14\ \ 
For the dispersion relation of a scalar in the near flat space limit, we need to evaluate the two sums 



XI (0 



n>0 
1/4 00 

9 V-/ -,Nn„, ^(4n+3){ng+0 



(D.22) 



n>0 
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Keeping the leading contribution to each of the coefficients 7714^+1, 'm.4n+3) we find from ()D.21|) 



yLO..^ ^ 1 V(_n»r(^)(4ra + 3)-^/^ r(n+|) (4„+3)g 



(D.23) 



(2^)'/\i^o r(n + i)r(f) 

To perform the two sums above, we introduce the integral representations 

r(|)(4n + l)-3/4 = / dss-V4e-(4"+i)^^ r(|)(4n + 3)-^/^ = / s^/^ e-(4"+3)% (D.24) 
Jo Jo 

and, after summing under the integrals, we obtain 
The asymptotics at ^ ~ —00 are given by 



For ^ ~ 00, the integrals (|D.25P are dominated by s < ^, leading to 

, 16 



(D.26) 



(D.27) 
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